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In recent years the study n f 

1 positron annihilation in 

solids has Become an important tn 

°I to study different physical 
properties of the solids. In . 

Xl cuiar the measurement of 

the angular correlation of posits 

°n annihilation radiation 
(ACPAR) has provided us with i m -nn>,+. 

v r tant data ah out the Permi 
surface (PS) and the electron mom ^ 

en tum distributions in metals 

and alloys. Recently different m 

a ohines have been developed 
to measure two-dimensional (2D) a n-p 

aoj? AR curves and such measure- 
ments in metals like Al, Cu an ^ 

^loys of Cu like Cu-Zn, Cu-Ni 
have shown that 2D-ACPAR technirm 

^ e is a powerful method to 

study the PS parameters. The 

on Cu, in particular, has 

shown that it is essential to oerf 

y • CI °rm a band structure calcula- 
tion of the two-photon momentum d-; 

uistributions in a transition 

metal for computing the SD-^qp*^ 

curves. By comparing the 

theoretical and experimental 2D-inr, 

"AUPAR curves one can obtain 

considerable details about the h 

a3a d structure and PS of 
transition metals. In the c&sp ^ 

01 a 3d transition metal like 
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vanadium (?) which, has a Fermi level lying among the d-bands, 
the hand structure and the FS are sensitive to parameters like 
crystal potential, exchange effects, self-consistency, electron 
configuration etc. used in the hand structure calculation. One- 
dimensional ACPAR measurements in the single crystal of V have 
already shown that the hand structure and FS characteristics of 
V give rise to fine-structure in these curves. It is, therefore 
expected that a 2D-ACPAR measurement in V would provide us with 
a greater wealth of data which can he used to test hand 
structure calculation. Keeping this aim in mind we have 
calculated the hand structure, FS, two-photon momentum distri- 
bution and the 2P-ACPAR curves for four sets of single crystal 
directions in V. The calculated 2D-ACPAR curves show fine 
structure which is analysed in terms of the hand structure and 
Fermi surface. The results are extended to predict the 2D-ACPAR 
curves in vanadium deuteriae ( VDq r^) within the framework of 
rigid hand protonic model. 

In Chapter 1 a general introduction to the subject of 
positron annihilation in metals is presented. After describing , 
the general characteristics of the phenomenon of positron 
annihilation the experimental techniques used to study positron 
annihilation in metals are described. Since the present work 
deals with the study by 2D-ACPAR curves, various experimental : 
arrangements used to measure the 2D-ACPAR curves are described 
in some detail. An outline of the theoretical calculation of 
positron annihilation in metals is given and the methods used | 


for comparing theory and experiment are discussed. At the end 
of Chapter 1 we have briefly discussed the motivation behind 
the present work. 

A review of the theoretical and experimental studies of 
the band structure of V is given in Chapter 2. The experimental 
study of the FS and momentum distribution in V made by different 
workers using the techniques of Fermiology , Compton scattering 
and positron annihilation is also reviewed. 

The method for the band structure calculation used in 
the present work was Hubbard's fast approximation scheme for 
transition metals. In Chapter 3 we have outlined various steps 
involved in this method and have discussed how Hubbard's fast 
approximation scheme has been modified by Mijnarends for the 
calculation of the electron wavefunctions and momentum distri- 
butions. Different methods used for the calculation of the 
positron wavefunction in metals are described and the correc- 
tions to the momentum distribution due to the many-body effects 
are briefly discussed. The Chapter concludes with an outline of 
the calculation of two-photon momentum distribution in metals. 

In Chapter 4 we have presented the results of our calcula- 
tions for the band structure* FS* two— photon momentum distribu- 
tion and 2D-ACPAR curves in V. Numerical details of the 
parameters used in our calculation are given. The results are i 
presented in terms of various tables and figures and are 
compared with the results of other calculations wherever 
possible. The calculated two— photon momentum distribution and 



the 2D-ACPAR curves are analysed in terms of the hand structure 
to show how different sheets of the FS of V can he studied with 
the 2D-ACPAR techniques. Although we have calculated the 2D- 
ACPAR curves, N(p ,p z ), for the four crystal directions A : p z = 
[OOl], p x = [100], B : p z = [001], p x = [110], C : p z = [Il2], 
p x = [ill], D : p z = [lio], p x = [lio] all in the Py = 0 plane, 
we have presented pictures of the PS topology in some of the 
p 7^= 0 planes parallel to the (100 ) and (110) planes to point 

‘ V 

out the importance of the 2D-ACPAR curves for p ^ 0. Our 
hand structure (2D-ACPAR) results are compared with those 
calculated hy a simpler modified Wigner-Seitz model to 
illustrate the point that hand structure calculations alone can 
explain the fine structure in the 2D-ACPAR curves in a transi- 
tion metal like V. 

Recently the electronic properties of metal-hydrogen 
systems have invoked considerable interest. In particular the 
hydrides and deuterides of the group V B metals V, Nh and Ta 
have been studied extensively, both theoretically and experi- 
mentally. In Chapter 5 we have extended the results cf our 
calculations for pure V (Chapter 4) to vanadium deuteride 
(YDq nsing the rigid hand protonic model. The results for 
VDq ^2 are discussed in terms of the changes brought in the 
hand structure and FS of V hy the introduction of deuterium. 
These results show how the 2D-ACPAR method can he very useful 
to study the hand structure and PS of VD . 
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The results obtained in the present work provide an 
understanding to the band structure effects in a 2D-ACPAR 
curves in V. In view of the similarity in the band structure 
and FS of the i soelectroni c metals V, Fb and la, the present 
results can be extended to understand 2D-ACPAR data in 1Tb and 
Ta. It is hoped that these results have emphasized the 
usefulness of the emerging 2D-ACPAR method for the study of 
transition metals and have provided guidance for designing 
(i.e. which ^p x ,Py,p z ^ directions to choose) a 2D-ACPAR 
experiment with V and VD . At present there is no 2D-ACPAR 
data reported for V and VD . The present results can provide 

” -A- 

a future basis for the comparison between experiment and 
theory. This work has also provided an opportunity to the 
author to learn about the interpretation of the 2D-ACPAR data. 




Chapter 1 


INTRODUCTION 


1.1 Introduction 

The determination of electron wavefunction is essential 
in understanding the electronic properties of atoms, molecules 
and solids. In addition to the other ground state properties 
of solids, the determination of the charge, spin and momentum 
distribution of electrons has recently become an interesting 
subject for experimental and theoretical research [l,2]. In 
particular the band structure and Fermi surface (PS) of metals 
is being increasingly investigated through the determination 
of the electron momentum distribution (EMI) of the metal. 

Three microprobes photons (Compton scattering) [ 3 ], positrons 
(positron annihilation studies) [4,5] and electrons (electron 
Compton scattering, (e,2e) scattering) [3,6] can be used to 
determine the EMD in metals [l,2]. 

The EMD, P(p), is defined as 

p(p) = lx(p)l 2 ( 1 . 1 ) 

where x(p) is the electron wavefunction in the momentum 
space and it can be obtained from the wavefunction, 4(r), in 
the position space through the Dirac transformation 

x (p) = / exp(-ip.r) dr 


( 1 . 2 ) 
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For a meaningful comparison between theory and experiment 
it is necessary to have a band structure calculation of the 
EMD in metals. Such a comparison can allow us to examine the 
effects of band structure on the EMD and to test the validity 
of the results of a particular band structure calculation. 

Among the different experimental techniques i.e. Compton 
scattering, positron annihilation and electron scattering, the 
Measurements employing two-dimensional (2D) angular correlation 
of positron annihilation radiation (ACPAR) are very useful in 
studying the F3 of metals. The present thesis describes the 
results of our calculations of the momentum distributions in V 
as studied by 2D-ACPAR methods. The band structure and FS 
topology for V is one of the most interesting among the 34 
transition metals. The results discussed in the present 
thesis provide suggestions for designing a 2D-ACPAR experiment 
to study the band structure and FS of V. These results are 
then extended to examine the YD system which has assumed 
considerable importance in the last few years. 

The plan of the present Chapter is as follows. The 
characteristics of positron annihilation process are discussed 
in Sec. 1.2 while the experimental techniques involved in the 
study of positron annihilation are described in Sec. 1.3* 

The theory for the angular correlation of positron annihila- 
tion radiation is outlined in Sec. 1.4. In Sec. 1.5 the ; 

methods used for comparing theory with experiment are described; 
and in Sec. 1.6 the motivation of the present work is discussed; 
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1 . 2 Characteristics of Positron Annihilation 

In recent j r ears diverse problems in the physics of 
condensed matter are being studied by using positron annihila- 
tion techniques [-4,5]. In these experiments low energy 
(0-700 keV) positrons from laboratory sources like radioactive 
itfa, ^Cu or- ? Go are allowed to annihilate in the sample 
under study and observations are made on the timing or 
energy or direction characteristics of the annihilation 
gamma rays. These low energy positrons, upon entering the 
material medium, slow down to thermal energies losing their 
energy through ionization and excitation collisions. The 

time taken for this thermalization of positrons is very short 
—12 

( % 10 s or so) compared to the lifetime for their annihila- 
tion (x > 10 ^ s). These features are shorn schematically 
in Tig. 1.1. 

A thermal! zed positron undergoes an annihilation with 

an electron by emitting a number of photons which will carry 

2 

an energy equal to the combined rest mass energy 2m Q c 
( Ail. 02 MeV). The sum of the kinetic energies of electron (few 
eV) and positron ( ^0.02 eV) is very low compared to the 
above energy (1.0 2 MeV). The probability for the various 
modes of the annihilation of the thermalized positron-electron 
pair is determined by various selection rules, conservation 
laws and other considerations using the invariance properties 
of quantum electrodynamics [7]. figure 1.2 shows the 
Teynmann diagrams for the most probable annihilation process 




Feynmarm diagrams describing one 
three-photon annihilation modes. 
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i.e. one, two and three photon emission. The ratio of the 
cross section for the one, two or three photon annihilation 
process is given by [ 4 ] * a and CT ( ly )A( 2 Y ) " 

where a = 1/137, the fine structure constant. Thus the 
probability of the two-photon annihilation is much greater 
than that for one and three photons. 

One can write the expression for the cross section for 
the two-photon annihilation at low velocities, v, of the 
positron as [ 4 ] 

°(2y) = " r o C / T U.3) 

2 2 

where r Q = e /m Q c = classical electron radius. The equation 
for the probability (r_ ) of two-photon annihilation becomes 

^ I 

r 2y = °(2y) ™e = 11 r o c n e U.4a) 

where n g is the electron density at the positron site. The 
positron lifetime v (t = l/ r 2y) ^ e i n S independent of positron 
velocity can provide the information about the electron 
density (n g ) at the site of annihilation in the sample to be 
studied. 

Equation (1.3) for the cross-section must be 
modified in such a way that it includes the Coulomb inter- , 

action between the annihilating positron-electron pair because 
the above expression is the approximate form of the relativi- 
stic expression obtained by Dirac [8]. The actual electron 
density seen by a positron is not only at positron site but is : 
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also spread in the region in which Coulomb attraction is 
effective. Taking these considerations and the velocity- 
dependent electron density into account one can write the 
two-photon annihilation probability which depends on the 
positron velocity as given by [9] 

t, f \ 2 / m N 2 2ua/|3 

W v) = nr Q fa) i _~ e - ^-p 2W|3) 


( 1 . 


where jB = v/c. 

The above expression for the probability for the two- 
photon annihilation is no longer valid in real metals because 
the following two corrections will have to be made: (i) the 
screening of the charge of positron in the many-el ectron 
system and (ii) the Pauli exclusion principle. A more 
realistic treatment which includes these effects is given by 
Kahana [lO] and it will be referred to later on. Since the 
scope of the present work is limited to the two-photon annihi- 
lation in metals, we shall not discuss the phenomena of 
positronium formation etc. 

In metals the kinetic energy of the electrons is of the 

order of few eV while the positrons being thermalized ( ^0.02 

eV) can be assumed to be at rest. The kinetic energy and the 

momentum of the annihilating positron-electron pair is, 

therefore, mainly due to the electron. In the centre-of- 

mass system the total energy available to the two photons is 
2 

2m Q c (2ll*02 MeV) and they are emitted in exactly opposite 
(0 = 180°) direction each having an energy of 511 keV. The 


4b) 
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conservation of energy and linear momentum (p) of the positron- 
electron pair in the laboratory system produces two measurable 
effects. Firstly there is a deviation from collinearity in 
the direction of emission of the two photons and secondly 
there is Doppler shift in the energy of the annihilation 
photons. These effects are illustrated in Fig. 1.3 which 
shows how the deviation © from collinearity arises from the 
transverse component p of the momentum p. The angular devia- 
tion 9 is related to p^ by the relation 

© - Pj/ m 0 c (1*5) 

and for the typical values of the electron momenta in metals 
© is of the order of few milliradians. The longitudinal 
component p. . of p causes the energy of the photon to shift 
from its value E = m Q c (i.e. Doppler broadening). It is 
seen from Fig. 1.3 that this energy shift AE is given by 


AE = c p /2 


( 1 . 6 ) 


We thus see that the measurement of the angle © and/or 
energies of the annihilation photons can provide useful infor- 
mation about the electron momentum distribution in metals. 

The details of the experimental methods using the above princi- 
ples are described in the next section. 


1 .3 Experimental Techniques 

1.3*1 Measurement of positron lifetime spectra and Doppler 
broadening of annihilation radiation lineshape 

We have seen above that the lifetime of positron is an 




8 



Fig. 1.3: 


The vector diagram of the conservation of momentum 
of the annihilation pair { p ) in the two«pho ton 
annihilation. The transverse and longitudinal 
components are denoted by J- and |i respectively. 
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important observable which provides interesting information 

about the sample under study. This is the most frequently 

used method in positron annihilation studies and it consists 
22 

of using a Na source sandwiched between the samples. The 

2 2 

1280 - keV gamma ray emitted by Na is almost simultaneous 
with the emission of a positron and it signals the ’birth' of 
positron. The 511 - keV gamma ray accompanying the annihila- 
tion of positron signals the 'death' of the positron. The 
lifetime measurement then consists of using the 1280-keV 
gamma ray as the ’start' signal and the 511-keV gamma ray as 
the ’stop’ signal in a time-to-amplitude converter system 
as in a typical nuclear lifetime spectrometer and of recording 
the time spectrum. Details of the experimental arrangements 
involved in a lifetime spectrometer are described in the 
literature [5,llj. 

Measurement of the Doppler broadening of the annihila- 
tion radiation lineshape is made with a high resolution 
gamma ray energy spectrometer involving a HP G-e or G-e(li) 
detector, research amplifier and a data acquisition system 
(e.g. multichannel analyser), as described elsewhere [2,5,ll]. 
There are also some multiparameter methods which determine 
more than one parameter and these are described in the litera- 
ture [5,ll]. Since the present work does not deal with the 
above methods we have not discussed them in detail. 
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1.3*2 Measurement o f the angular correlation of -positron 

annihilatio n radiation 

We have seen earlier that the measurement of the 
angular correlation of positron annihilation radiation (ACPAR) 
leads to the determination of the momentum distribution of 
the annihilating positron-electron pair. The geometry of a 
typical ACPAR measurement is shown schematically in Pig. 1.4 
where the double arrow represents the electron momentum vector 
p. It was shown in Sec. 1.2 that the transverse component 
p^ = (Py>P z ) i s responsible for the deviation from collinearity 
of the two photons. If one of the photons is detected by a 
detector (Pig. 1.4) the second time-coincident photon will 
strike the other detector plane C2 at a point whose coordinates 
are y 1 9 y 1 and z L (since both 9 and © z are of the 

order of few milliradians) where 1 is the distance between the 
annihilation centre 0 and each of the counter planes. The 
two angles 9 and © are related to the momentum components 
Py and p 2 through the relations © z = V z / m 0 ° and © y = Py/m Q c. 

The longitudinal component, p , gives rise to a Doppler energy 

-A. 

shift &E which has a typical value of AE = 1-3 keV away 
from Ey = m 0 c or ^ 511 keV. In order to achieve high 
counting rates the gamma ray detectors in the ACPAR experiment 
have high detection efficiency but poor ( ^10 '/, at 511 keV) 
energy resolution. As a result the Doppler energy shifts AE 
are not resolved and the momentum component p x is undetermined 
or is integrated out. There are two main variations in the 
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geometry of the ACPAR experiment. In the so-called ,r long-slit 
(IS)’’ (or wide slit) geometry the detector on each side has an 
extended length along the p y -axis in order to enhance the 
coincident counting rates. In such a measurement the informa- 
tion about the p -component is also integrated out and the LS 
ACPAR measurement yields the coincidence rate proportional to 

CD 

N (p 2 ) = / /“p(p)a p x dp y ( 1 . 7 ) 

— 05 — CD 

where P(p) is the two-photon momentum density (TPMD) and 
where we have assumed that the horizontal angle subtended by 
the detectors is much larger compared to the angular resolution 
of the detectors along p z . 

By sacrificing the coincident counting rate one may 
restrict the dimensions of the detectors along the p -axis 
and employ the so- call ed* point slit (PS)" geometry . In this 
type of geometry one determines the p -component also and 

V 

the ACPAR measurement yields a two-dimensional distributions 

co 

N(p y >P z ) = / p(p) dp x (1.8) 

Having discussed the principles behind the IS and PS 
geometries we shall now describe the actual experimental arra- 
ngements used. 

A schematic arrangement of a typical ACPAR set-up ; 

employing the LS geometry is shown in Pig, 1.5(a). Positrons ; 
emitted from a radioactive source like (or ^Cu or ^^Co) ; 
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are focussed "by an electromagnet onto the sample under study. 

The annihilation of positrons in the sample results in two gamma 

rays (E^ = 511 keV) which follow collimated paths towards the two 

gamma ray detectors placed at about 3 meters on each side. The 

gamma ray detectors consist of Nal(Tl) scintillators cylindrical 

in shape with about 5 cm diameter and about 30-80 cm in length. 

The axis of each detector is horizontal and the entrance window 

is also horizontal and it subtends an angle of about 0.5 mrad 

at the sample. One of the detectors is fixed while the other 

is movable along the arc of the vertical circle centered at the 

sample. The measurement consists of recording the coincidence 

counting rate N(9 = p/m c) as a function of 9„. It is seen from 

eq. (1.7) and Pig. 1.6 that the LS-ACPAR measurement counts the 

number of electrons enclosed in a slice (parallel to the p p 

x y 

plane) in the momentum space when this slice is moved parallel 
to itself. In case of a free electron gas (i.e. Eermi sphere) 
the LS-ACPAR curve IT(9 = p/m c) should be described by an 
inverted parabola. Actual measurements in simple metals (e.g. 

Li whose LS-ACPAR curve is shown in Pig. 1.5(b)) indeed show 
that the ACPAR curve is an inverted parabola superimposed on a 
broad Gaussian-like distribution. The parabolic part is attri- 
buted to the nearly free (or conduction) electrons while the 
Gaussian part arises out of the annihilation with inner core 

t 

electrons. A sharp break at the Fermi momentum p = pj, is also 
seen. The nature of the LS-ACPAR curve in a transition metal ; 
is more complicated in that it is anisotropic and has a : 
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substantial contribution from high momentum components. 

The measurement of two-dimensional ACPAR curves, described 
by eq. (1.8), can be made by using different types of experimental 
arrangements. In view of the interesting information obtainable 
from the 2D-ACPAR measurements, much effort has been put in 
developing machines to measure the 2D-ACPAR curves [12,13]. Some 
of the machines developed so far are multiple Nal(Tl) detector 
system, multiwire proportional chamber system, Nal Anger camera 
system and multiple solid state detector system. We shall now 
briefly describe each of these arrangements. 

A) Multiple Nal(Tl) detector system 

Basically this is point-slit (PS) system whose detection 

efficiency is increased by using a two-dimensional array of 

Nal(Tl) detectors on each side of the source and specimen 

assembly (Pig. 1.7). Such a machine was first built by Berko 

and coworkers at the Brandeis University, U.S.A. and it employed 

22 Nal(Tl) detectors subtending an angle of 0.5 mrad x 2.2 mrad 

at the sample. A newer version of this machine has recently 

been built at the Brandeis University and it consists of two 

sets of 32 Nal(Tl) detectors on each side, with the plane of 

detectors being 10 meters away. Each of the detector is a 

Nal(Tl) scintillator having a length of 5 cm and diameter of 

3.8 cm coupled to RCA 6199 photomultiplier tube. Collimation : 

of the gamma rays is made through rectangular lead point slits 

2 

subtending an angle 0.5 x 1.5 mrad at the sample. The number j 





The results of 2D-ACPAR measurements for 
Cu reported by Berko [12] for the crystal 
directions p^ » [101 ] * py « [121 J * 
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of (Py»P z ) pairs obtained by 32 x 32 cross coincidences is 109 

and more intermediate points can be obtained by moving the 

counter assembly. Collection of data is done in steps of 0.2 

mrad along the p direction in the range p = - 40 to + 40 mrad 

and in steps of 0.3 mrad along the p direction in the range 

Py = 0 to 20 mrad. The source of positrons is the radioactive 
58 

isotope Co and with a 500 mCi of such a source the measured 
counting rate for Cu sample is 0.2 counts/sec per. counter pair 
at Py. = p z = 0 (with a total coincidence rate of 50 counts/sec). 
The counting period required to collect a full 2D-ACPAR curve 
with a 0.4 X accuracy is of the order of 10-14 days. The 
results of such 2D-ACPAR measurements for Cu reported by Berko 
[12] for the crystal directions p z = [lOl], p y = [121 ] are 
shown in Pig. .1.8 and they bring out the effects of the Fermi 
surface of Cu. 

B ) Multiwire proportional chamber system 

The design of this system is, in principle, similar to 
the apparatus built by Howell and Osmon [14] who measured coin- 
cidence counting rate between a spark chamber and a point-slit 
Hal(Tl) detector to determine the 2D-ACPAR curve. The useful- ! 
ness of the spark chamber machine was limited by the poor (^2 X ; 
detection efficiency. During the last few years the group at 

I 

CERN, Geneva and the University of Geneva, Geneva, Switzerland ! 

I 

have developed a 2D-ACPAR machine consisting of multi wire 
proportional chambers (MWPC) having high density converters to 
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increase the efficiency [l2j. In this macliine the annihilation 
photon entering the chamber is detected through the photo- 
electron formed by it after interaction with the high-Z 

2 

converter. The chamber has a sensitive area of 10 x 10 cm and 
the detection efficiency for the 511-keV gamma ray is about 12 % 
The two M'WPC detector systems are placed 10 meters apart and the 
data are recorded by a CAMAC system coupled to a PDP 11/40 
computer system. The total coincidence counting rate is 10 

coinc/sec while the effective angular resolution is 0.75 x 0.68 

2 22 
mrad when a 30 mCi radioactive source of Na and a copper 

p 

sample (10 x 1 nmr) was used. The 2D-ACPAR curve for Cu 
obtained with this MWPC system is similar [l 2] to that shown in 
Pig. 1.8. 

C) Nal Anger camera system 

The position sensitive or Anger camera was developed 
many years ago [l 5 ] and has been used for medical imaging for 
many years. In principle the Anger camera consists of a large 
Nal crystal optically coupled to several photomultiplier tubes 
and the localization of the gamma ray detection point (or 

I 

position) is made by analysing these signals. A 2D-ACPAR 
machine using a Nal Anger camera has been built by West and 
coworkers at the University of East Anglia, U.K. [12]. The 
Anger cameras used in this machine consist of two cylindrical 
Nal scintillators 50.8 cm in diameter and coupled optically to • 
a closely packed array of 37 photomultiplier tubes each having 
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a diameter of 7.8 cm. The source and sample assembly is at the 

centre of two 'cameras' which are operated in coincidence. The 

distance between the two cameras can be changed from 2 to 20 

meters. A typical total coincidence counting rate of the order 

of 120 coinc/sec (using a 20 mCi Na source and an aluminium 

sample in the range -22 to +22 mrad) and an overall angular 

2 

resolution of about 0.6 x 0.6 mrad is reported. 

D) Multiple solid state detector system 

It was pointed out by Triftshauser in 1971 that a 2D-ACPAR 
machine can be constructed by using a position sensitive multiple 
de-detector system. Such a 2D-ACPAR machine has been constru- 
cted at the Kernf or s chung sanl age , Julich, federal Republic of 
Germany [12]. This instrument consists of five point-slit Osl 
detectors on one side and a high purity germanium (Ge) position 
sensitive detector on the other side with a coincidence between 
the Csl detectors and the elements of the Ge detector. The 
dimensions of the Ge detector are 90 mm x 60 mm x 10 mm with a 
checkerboard type of design Cl 2 ] . In a similar design, the 
research group of the Institute of Physical and Chemical Res ear cl 
Saitama, Japan has developed a 2D-ACPAR machine using large ares 
planer Si(Li) detectors [12 ] . The data by these two types of 
2D-ACPAR machine have not yet been reported in detail. j 

In the above we have described, in brief, the significant 
features of various experimental arrangements employed to 
measure the 2D-ACPAR curves N(p„,p„). In the next section we ; 

J & 
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shall outline how the N(p.p_) curves can be calculated 

y z 

theoretically. 

Before we conclude this section we wish to point out 
some distinct advantages offered by the 2D-A0PAR technique. 

In a pure metal the study of Fermi surfaces can be carried out 
by other conventional methods such as the de Haas -van Alphen 
effect, the magnetoresistance effect etc. with a precision 
higher than that offered by positron annihilation techniques. 
However these conventional methods cannot be used for the study 
of nondilute disordered alleys and high temperature phases of 
ordered materials because of the short mean free path of 
electrons under these conditions. Although there are other 
methods like Compton scattering, X-ray diffuse scattering, 
inelastic neutron scattering, Faraday effect etc. available for 
studying samples with solute concentration *1 /. , the superior 
momentum resolution and the possibility to measure 2D-ACPAR 
curves make positron annihilation a powerful tool for studying 
nondilute alloys. Many such interesting studies have been 
reported recently [5,12]. 

The positron method has a serious drawback in that the 
sample under study should be free from defects. This restri- 
ction arises from the fact that a defect present in a metal 
can 'trap' the positron changing the TPMD significantly [ 5 ]. j 
In our discussions throughout the present thesis we have limited;' 

I 

ourselves to pure metals without any defects, i 
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1 . 4 Theory of Angular Correlation of Po si tron Annihilation 
Radiation ; Independent Particle Model 

It was seen in the previous section that theoretical 
determination of “both one- dimensional ACPAR curve, N(p ) by 
eq. (1.7) and the two-dimensional ACPAR curve, N(py,p z ) by 
eq. (1.8) require calculation of the two-photon momentum 
density (TPMD) P(p). Such calculations of the TPMD, P(p), 
are done within the framework of the independent particle 
model using a standard band structure method [4,5]. The TPMD, 
p(p), is the probability that the photon pair carries away 
a total momentum p and in the independent particle model one 

can write P (p) in a metal as 

occ 

P(p) = Const £ f(k,j) | A^(k,p) | 2 6(p-l-K i ) (1.9) 

3,k,i 

where f(k,j) is the occupation number such that f(k,j) = 1 

-y 

for an occupied state and f(k,j) = 0 for an unoccupied state 

-y -y , 

of the electron at T = 0 K, the overlap matrix element A-(k,p) 

J 

is written as 


A i (2,p) = / exp (-ip.r) 

J Cell 

where .(?) is the electron 

k, 0 

and band index j and 4 (r) is 

-y 

the state k, = 0 in its lowest 

-r 

reciprocal lattice vectors. 


4 + (r) 4 (?) dr (1.10) 

k, 3 + 


wavefunction for wavevector k 


the positron wavefunction for 

~y 

energy state and K^'s are the 
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It was pointed out by Berko and Plaskett [ 16 ] that 

-V 

because the crystal momentum k differs from the real momentum 
p by , an electron with a wave vector k will contribute to 
the TPMD not only at p = k but also at p = k + 2^ (Fig. 1.9) 
with amplitudes A^(k,p) through eq. (1.9). The contribution 
at p = ic (with = 0 ) are termed 'normal' contributions while 
the contributions at the high momenta p = k + (with ^ 0 ) 
are termed as Umklapp contributions or the high momentum 
components (HMC) . The summation in (1.9) is carried over the 
occupied states and as a result the TPMD, p(p), is expected 
to show discontinuity not only at p = icj, (where k-p is the 
Fermi momentum) in the first zone but also at the higher 

-f ->• ->■ 

momenta p = kj, + Besides these discontinuties which 

arise out of the Fermi surface, the general shape of the 

TPMD, P(p) , is determined by the electron wavefunction 

.(r) via eq. (1.10). In the case of a transition metal 

like V which has partially occupied 3d bands, the TPMD P(p) 

becomes anisotropic because of this dependence on $•+ .(r). 

ic, 3 

There is also a symmetry selection rule [ 17 ] which allows only 
bands of certain symmetry to contribute to the TPMD and this 
aspect will be discussed in detail in Chapter 4. 

In order to illustrate the important features of the 
TPMD, p(p) , in real metals we take the example of one- 
dimensional nearly-free-el ectron (NFE) and tight binding (TB) 
model for metals. Three cases corresponding to different 
locations of the FS with respect to the Brillouin zone (BZ) 
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boundary are shown in Fig. 1.10. In the case (a) the FS is 
within the first BZ while in the case (c ) the FS is in the 
second BZ (using the extended zone scheme). For the case (b) 
it is assumed that the FS coincides with the BZ boundary. The 
schematic trends in Fig. 1.10 show how the FS breaks appear, 
for the case (a) in the p-distributions at higher momenta. Cae 
(b) illustrates how the ^-distribution becomes continuous as 
soon as the n- distributions touch the BZ. In the case (c) one 
notices that the FS discontinuity is more sharp in the second 
BZ than in the first BZ. 

Coming back to eqs. (1.9) and (1.10) we see that the 

calculation of the TPMD, p($ 0, should be preceded by a band • 

structure calculation of the metal so that the FS as well as 

the electron wavefunctions ifu .(r) are determined. This 

k , 0 

stage of computation involves the standard problem of band 
structure calculation and it is discussed later in the thesis. 
Next problem is the determination of the positron wavefunction 
in the metal and this is treated in Sec. 3.4. Having calcula- 
ted the electron and positron wavefunctions one can determine 
the TPMD by the help of eqs. (1.9) and (1.10). Finally the 
calculation of the ACPAR curves can be performed with the help 
of eqs. (1.7) and (1.8). In particular the calculation of the 
2D-ACPAR curves N(py,p z ) require determination of theoretical 
TPMD, p(p), at several p-points before one can evaluate 

co 

N (Py» P z ) = / P'(f ) dp x ( 1 . 8 ), 

00 I 
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through a line integration. 

It should be pointed out that the positron-electron 
many body effects affect the positron lifetimes very seriously 
but the ACPAR curves only marginally. These corrections will 
be discussed in Sec. 3.5. In the next section we shall out- 
line how one can compare the theoretical and experimental 
results for the ACPAR measurements. 

1 . 5 Comparison Between Theory and Experim en t: Reconstruction of 
Momentum Distributions . 

In the last two sections we have seen how the ACPAR 
curves can be measured experimentally and how they can be 
calculated theoretically. There are two main approaches for 
comparing the experiment with theory. In the first approach 
one calculates the theoretical TPMD, P(p) , by using a band 
structure method and then uses eqs. (1.7) or (1.8) to obtain 
the ACPAR curves. Por an accurate comparison one will have to 
apply correction for the finite angular resolution of the 
apparatus as well as for the finite length of the long slit 
in the LS-ACPAR apparatus [ 5 ]. Such a calculation of the ACPAR 
curve from the TPMD, p (|T) , is rather straightforward for the 
2D-ACPAR curves because it involves a line integral (1.8) 
which can be computed by determining p(p) on a rectangular 
momentum grid in the (p p ) plane p = const. This procedure, 
followed by us in the present work, will be illustrated in 
Chapter 4. As for the 1D-ACPAR curves are concerned the 

surface integral involved in (1.7) can be calculated either : 
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by a direct integration in the momentum space [18] or by the 
help of Fourier -Hank el transforms [l9,20]. 

In the second approach one attempts to reconstruct 
the TPMD from ACPAR curves measured for several single crystal 
directions and compares the reconstructed TPMD with the 
theoretical results for p(p). We shall now briefly describe 
the methods used for the reconstruction of the momentum 
distributions. 

The method most commonly used for the LS- ACPAR curves 
is due to Mijnarends [20] and is based on the Fourier-Hankel 
transforms. In this method the momentum density P(p) is 
expanded as 

P(?) = l P„ (p) (1.11) 

where describes the orientation of p while Av ( are 

the lattice harmonics which form an orthonormal set of invariant 
linear combinations of spherical harmonics of order a. The 
label v is used to distinguish different harmonics of the same 
order and symmetry. Only the harmonics having z= 0,4,6,8,10, 

12 (twice), 14 will contribute for cubic lattices like V 

(bcc). Two orthogonal co-ordinate systems R r (fixed to the 
crystal) and R (fixed to the apparatus) are used. The 
lattice harmonics are defined in the system R = ( £ , t) , c ) 
while the z-axis is selected in the system R = (x, y, z). 

By a suitable rotation through the Euler angles (a, p, y) 

the systems R~ and R^ can be made to coincide. j 
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The lattice harmonics F £v ( fi^) can be further expanded 
in terms of the spherical harmonics Y £m ( by the following 

relation 

^ = ^ a £mv Y £m^ n p^ (1*12) 

x m ^ 

Following Hijnarends the 1S-ACPAR curve Np a (p z ) is expanded as 

= 2E L Fj ' u (P,Il) s l» iv z' > (1 - 13) 

where the angles (|3,a) denote a particular orientation of 
the single costal specimen and 

K> 

g £v (p ) = f p£v(p) \ (p z /p) P d P (1.14) 

p z 

where P £ are the Legendre polynomials. 

Equations (l.ll) to (1.14) constitute the scheme 
for the reconstruction of p(p) which proceeds as follows. A 
number of LS-ACPAR curves, N(p ), are measured for different 
single crystal orientations (oc, (3). The expansion (1.13) is 
truncated after a finite number of terms, m, such that m<_n 
where n is the number of orientations (a,|3) for which measure- 
ments are made. By fitting experimental data, N(p ), in eq. 
( 1 . 13 ), one solves for g„ v (p z ) and then obtains p £v (p) from 
the Fourier-Hanhel transform f 

P j 

= -“ / g” (2) £,(z/p) d 2 ( 1 . 15 ) 


P£v(p) 


o 
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1 

P 


r d §£v( p) £(i+l) 

1 3F 2p 


g £v (P) 


P 

+ ~o / S (z) P£' (“) dz] (1.16) 

o ^ P 

Jhe TPMD, p(p) , is then synthesized from the p^ v (p) using 
(l.ll). This method has "been used for the reconstruction of 
the EMD in several metals [3,5]. One problem inherent in 
this method is the arbitrary choice of m while truncating 
(1.13). 

The problem which is more closely related to the 
subject of the present thesis is the reconstruction of the 
TPMD, p(]?), from the 2D-ACPAR curves, R(p„,p_), measured 

j * 

experimentally for a set of {p x , Py, p z } directions. There 
are no reported measurements on the 2D-ACPAR curves for V in 
the literature. Nevertheless, the theoretical TPMD curves, 
p(p), calculated by us and presented in Chapter 4 can provide 
a theoretical basis to compare the reconstructed TPMD once 
experimental data become available. 

There are several methods suggested for reconstructing 
TPMD from the 2D-ACPAR curves [ 21-23 J. Majumdar [ 21 ] was 
the first to suggest a method which is based on Mellin 
transforms and which was aimed directly on the problem under 
consideration. 

Let us consider, for simplicity, the measurements in 
the 1 ' fixed configuration” such that the axes of cubic crystal 
coincide with the axes of measuring system. In the co-ordinates 
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(R ) of the measuring apparatus we can then write p(j$) as 


p V 

p(j7) = l -t Zm Pa (cose) cos m 0 P £ v(p); ( m f.0 ) 


ilvm 


••• (1.17) 


One can write eq. (1.8) in the polar co-ordinates (R, 0) as 


N(p y ,p a ) = g(R, 0) = l g m (R) cos in 0 


m 


(1.18) 


where g (R) is the Fourier cosine transform of the point-slit 
geometry data given by 


g m (R) 


\ v Urn C( R ) 


(1.19] 


with 


(R) = / [(i - R^P *)*] p Av (p) (i - R 2 p 2 )~* d P 


,2 - 2 . 


R 


( 1 . 20 ) 


,m 


To find out p£ V (p) one has to solve eq. (1.19) for f £v 


and 


then eq. (1.20). for p p (p). Magumdar suggested constructing 
the Mellin transform (Re s >1) 


?m 


(s ) = / R S “ 2 f“ u (R) dR 


( 1 . 21 ) 


for writing the solution of eq. (l.20) as 


£v 


{ P> = Hi 7 


cr+i ” 2 1+m r(l+m) r (s+m)) 

__ - 

o'-i«° P 


Pjfv (s) 


Tt 2 r(w (s+m+l)) 


X 


3P 2 (£+m+l), - Tj(il-m), ^ (s+m+l); 1+m, ^ (s+m); l) 


( 1 . 22 ): 



31 


where o > 1 . 

m 

In eq. (1.19) fj, v 's are again underdetermined as "before 
but it is easy to find out the maximum number of ^-values to 
be used because the m-values in using eq. (1.18) with the 
experimental data are known. It may be noted here that any 
information about the m-values cannot be obtained using the LS- 
ACPAR data. Above procedure has been discussed in detail by 
Majumdar [2l]. His method is more systematic and gives a 
unique criterion for deciding where to truncate the expansion 
(1.17). 

Cormack [22] has suggested a method to find a real 
function in a finite region of a plane if the line integrals 
along all the straight lines intersecting the region are given. 
The solution given is applicable to three problems of interest 
in radiology and radiotherapy (i) the determination of a 
variable X-ray absorption coefficient in two dimensions j 
(ii) the determination of the positron annihilations when 
there is an inhomogeneous distribution of positron emitter in 
the matter • and (iii) the determination of variable density 
of matter with constant chemical composition using energy loss 
of charged particles in matter. In particular these methods 
can be used for the reconstruction of the TPMD from 2D-ACPAR 
measurements. 

Recently Pecora and Ehrlich [23] have reported a method 
to reconstruct the momentum distribution from point- si it 
geometry data, particularly the crossed slit geometry. They 
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have used the following expression 

N a (P z ) = / p(p x > P y = °> P z ) d P x (1.23) 

— co 

for the experimental data in a plane with the z-axis in the 
co-ordinate system of experimental arrangement making an angle 
a with the reference direction in the plane. They have 
expanded the p(^) in a Pourier series in polar co-ordinates 
in the plane using the relation 

p(p) = l P m (p) (1.24) 

m=- 00 

where p s J]T| and 0 is the angle of p with respect to the 
reference direction. The p(jf) must be real and this results 
in the following condition 


R e[p m (p)3 = Pm(p) = P_m(p) 

and 

Im [ =- p _m(p) 


It is therefore required to solve for the coefficients P m (p) 

only for m >0. To find out IT (p„) one uses eq. (1.24) with 

OC * 47 z 

, D 

0 = a.+ 0 where a i is the angle between the z-axis and the 
reference direction while © is the angle between p and the 
z-axis. One can write 

«> ima. im0 

p(p) = I e J P m (p) e (1.25) 

-*■ co 

and 1 

1:0 ima. « ; 

N a ( ?z } = ^ e 3 / p m (p) cos m© dp x (1.26; 

U *4 00 l 


— oo 
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The eq. (1.26) can he rewritten as 


N a .(p 2 ) = l g m (p z ) e 

D m=- «= 


ima . 
tJ 


with 


g m (p z ) = 2 / 


■ / - ItI" Vp> T m<P z /P) 


P, 


( P ‘ -p;)- 


(1.27) 


(1.28; 


where the Chebychev polynomial of the first kind, T , is 
defined by 


T (x) = cos [m cos ^ x] 

Now one has to solve eq. (1.27) for g m (p z ) and then 
eq. (1.28) for P (p). The details about solving these 
equations are given by Pecora and Ehrlich [23]. They have 
tested this method for the hypothetical free electron gas 
and analysed the effect of errors in the data on the recons- 
tructed, p(p). They have also applied this method to real 
data for Cu and. tried to find out the information about 
various quantities like p-p. They found their results to be 
in remarkably good agreement with the theoretically predicted 
and reconstructed p(p) for Cu in the (100 ) plane. They were 
able to show that the data for a minimum number of crystal ■ 
orientations can provide rich information about, P(p). This 
method appears to be promising and future applications of 
this method for the reconstruction of the TPMP from experimental 
data should provide considerable interest. 
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1 . 6 Motivation Behind the Present Work 

In the previous sections we have pointed out that the 
2D-ACPAR method is well-suited for testing the band structure 
and Fermi surface of metals and alloys via the two-photon 
momentum distributions. Earlier work on the 2D-ACPAR results 
for A1 [24] and Cu [25,26] have shown that for a successful 
interpretation of the experimental data it is essential to 
perform theoretical calculations made by band structure 
methods. Such calculations are particularly important for 
those transition metals for which the 3d band is partially 
occupied and for which the d wavefunctions play a major role 
in building the TPMD through eq. (1,10). One such metal is 
V belonging to the group V B metals (i.e. 7, Nb and la). The 
band structure and PS of these three metals 7 , Nb and Ta are 
similar and they depend on the crystal potential, exchange 
character, self-consistency, electron configuration and the 
band structure method used in the calculations. The calculated 
as well the measured 1D-ACPAR curves in 7 have already shown 
a highly anisotropic nature which could be understood in 
terms of the band structure and Permi surface of 7. 

We therefore thought it interesting to calculate the 
TPMD and the 2D-A0PAR curves in V using a band structure 
method based on Hubbard’s approximate scheme. The choice of 
the crystal potential for metallic 7 for these calculations 
was based on the following considerations. Under the 
auspices of the International Union of Crystallography a 
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Special Commission on Electron Charge, Spin and Momentum 
Density was formed in 1972, In 1975 this Commission piloted 
a special project to determine the accuracy of the theoretical 
Compton profiles and X-ray scattering factors for metallic 
V, Our research group took part in this International 
project and has calculated the electron momentum distribution 
(for Compton scattering) and the Compton profiles [27] using 
Hubbard' s fast approximate scheme. These calculations used 
the standard lattice potential for metallic V specially 
constructed by Mijnarends [28] for this project. We have, 
therefore, used the same crystal potential for V and have 
extended the previous work on this project by calculating 
the two-photon momentum distribution (for positron annihila- 
tion) and the 2D-ACPAR curves in 7. 

The aim of the present work is to examine the effects 
of band structure on the TPMD and the 2D-ACPAR curves in V. 

Ac present there is let of interest in the electronic struc- 
ture of transition metal hydrides. The hydrides and deuterides 
of vanadium have been studied both theoretically and experi- 
mentally [29]. We have, therefore, extended our results for 
V to VDq ^ using the rigid band protonic model and have 

calculated the 2D-ACPAR curves in VD n nn . The results of 

u . 1 c 

the present work can give a better idea about which PS 

sheets can be investigated in a particular choice of 

{p x ,Py,P z > direction in a 2D-ACPAR experiment with V and 

VD . 
x 
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At present there is no available experimental data on 
the 2D-ACPAR measurements in V or VD . \ie hope that the 

-X. 

results presented in Chapters 4 and 5 can provide a theore- 
tical basis to understand future experimental results. 
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Chapter 2 


REVIEW OE THE THEORETICAL AND EXPERIMENTAL STUDIES 
OE THE BAND STRUCTURE OF VANADIUM 

2.1 Introduction 

The study of Eermi surface (ES) of metals and alloys 
plays an important part in the investigation of their elec- 
tronic properties. Determination of the ES topology is 
particularly interesting in transition metals having 
partially occupied d bands. In the past there has been 
considerable interest in the experimental as well as theore- 
tical study of the band structure and Eermi surface of the 
vanadium group of metals (i.e. the V B group is ©electronic 
metals, V, Nb and Ta) . The outermost d-shell (e.g. 3d for 
V) is partially occupied in all these three metals and the 
Eermi level, therefore, lies amongst the d bands thus 
resulting in a complicated topology of the ES. These 
interesting features of the ES topology can be determined 
by a comparison between theory and experiment. In case of 
the transition metals, the calculated band structure, Eermi 
surface and the electron wavefunctions are particularly 
sensitive to the band structure method, crystal potential, 
outer electron configuration, exchange parameter etc. used 
in solving the Hartree-Eock-Slater equations (in the 
independent particle model) 

= [ - 72 + E c (r) + V x (f)] *(r) 


H f(r) 


( 2 . 1 ) 
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where V (r) is the non-exchange potential experienced by 
the electron in the Coulomb field of all ions and all other 
electrons while V x (?) is the exchange potential which 
arises because the electrons have spin, are indistingui- 
shable and obey the Pauli exclusion principle. In eq. (2.1) 
we have used the atomic units : ii = 1, m = 1, e=l while 
the energies are expressed in Rydbergs. 

A detailed description of how eq. (2.1) can be solved 
by different band structure methods is given in several 
excellent reviews Cl— 5 ] and will not be attempted here. In 
the next section we shall briefly outline the different 
theoretical calculations of the band structure of metallic 
vanadium and summarise the experimental results for the PS 
topology of V. In Sec. 2,3 we shall summarise the experi- 
mental results on the electron momentum distributions in Y 
as investigated by the methods of positron annihilation and 
Compton scattering. 

2.2 Theoretical Calculations of the Band Structure of V 

The first band structure calculation of Y was 
reported by Mattheiss [6] who used a non-self-consistent 
(NSC) augment ed-pl an e- wave (APW) method and showed that 
the band structure of Y was affected significantly by 
the change in the relative population of the (3d) and 
(4s) shells. In a subsequent work Mattheiss [ 7 ] predicted 
the PS topology for V , Nb and Ta. His results showed that 
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the PS of vanadium group of metals consists of two hole 
sheets. The first one, arising out of the second band, forms 
a closed surface centered on the point r and is octahedral in 
shape. The second sheet, contributed by the third band, 
consists of two geometries : (i) jungle-gym (or JG-) struc- 
ture of hole arms along the <100 > directions and (ii) large 
ellipsoids around the points N. The PS topology deduced by 
Mattheiss [7] from the band structure of Fe, however, showed 
a neck along rN due to the third band, 

A band structure calculation by Snow and Waber [8] 
using a NSC APW method with full (a = l) Slater exchange 
favoured an outer electron configuration (3d)^ (4s)^ for V 
when a comparison with experiment was made. 

A NSC APW calculation of the band structure of V was 

T O 

made by Anderson et al. [9] who assumed a (4s; 

configuration and used two different values of the exchange 
parameters a = 1 and a = 2/3. These authors observed that 
the d-bandwidth decreased with increase in a. When their 
theoretical results were compared with the experimental data 
obtained by photo emission [lO] and de Haas-van Alphen effect 
[ll] studies the a = 1 type of exchange was preferred. The 
band structure of V was also calculated by Christensen for 
the (3d)^ (4^)^ configuration and these results agreed with 
those of Mattheiss [6]. Yasui et al. [12 ] made a self- 
consistent (SC) calculation on 7 using a modified tight 
binding method for the d electrons and an OPW method for 
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the s electrons. A comprehensive calculation of the hand 
structure of V was made by Papaconstantopoul os et al. [l 3 3 

rz p 

using a SC APW method and a (3d)^ (4s ) c configuration. These 
authors reported their calculations for two values of the 
lattice constant a Q : (i) normal lattice constant a Q with 
a = 1 and a = 2/3 and (ii) reduced lattice constant a = 0.95 a 0 
with a = 2/3. When the different theoretical calculations 
were compared to the experimental data obtained by photo- 
emission [lO] and de Haas van Alphen [ll] methods the SC 
calculation using a = 2/3 gave a good agreement. These 
authors also found that the band structure and FS topology of 
V showed significant changes on changing the lattice parameter 
from a Q to 0.95 a Q . 

Hattox et al. [14] have reported a SC APW band struc- 
ture calculation of V employing a statistical exchange X a 
[ 3 ] muffin- tin approximation. These authors have calculated 
cohesive energy, pressure and magnetization for several 
values of the lattice parameter. 

A detailed study of the band structure and FS of V 
has been carried out by Wakoh and Yamashita [15] using a SC 
KKR method. In order to match the experimental results on 
the dimensions and shape of FS these authors used a pheno- 
menological state-dependent potential having different 
energy dependence for the de and dy^ states. Their calcula- 
tions gave a good description of the data for FS topology 
and X-ray form factors. Using this state-dependent potential 
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these authors also calculated directional Compton profiles 
[16] and positron angular correlation curves [l7] and we 
shall discuss these results later on. 

The hand structure calculation of the group Y B metals 
V, Nb and Ta using SC APW method is reported by Boyer et al. 
[18] who found that the band structure of Y is very 
sensitive to the exchange approximation used. 

Recently the band structure of Y, Fb and Ta has been 
calculated by Alward et al. [l 9 3 using empirical pseudopoten- 
tial method. These authors compared the theoretical reflecti- 
vity data, density of states and Fermi surface topology with 
experimental results and found a general agreement. 

There are three band structure calculations of the 
electron momentum distributions (EMI*) in Y. Besides the 
work of Wakoh and Yamashita [16,17] mentioned earlier, 

Laurent et al. [20] have used a self-consistent LOGO method 
to calculate the Compton profile of V and other paramagnetic 
3d metals. A calculation by Kanhere and Singru [2l] gives 
the systematic s of the EMD in 3d metals including Y. 

Recently Wakoh and Kubo [22] have calculated the 
X-ray form factors by the APW method. They have estimated 
the charge asphericity in Y by using the ratios of the 
integrated intensities of X-ray scattering for reflection 
pairs. They find that the state-dependent potential constru- 
cted by them for Y (and which was fitted to the PS data) 
gives a closer agreement with the X-ray form factor data 
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than the state-independent potential [23]. 

It is clear from above that several authors have 
calculated band structure and other electronic properties 
of V. Although we have not discussed each of the se cal cul a- 
tionsin detail, we shall refer to their important conclusions 
in the next section as well as in Chapters 4 and 5. However, 
the results of the different band structure calculations 
reveal the following significant factors which determine the 
band structure of V : (i) The relative population of the (3d) 
and (4s) electrons assumed in the calculation, (ii) the 
exchange parameter, a, used, (iii) whether the calculation 
was performed in a non-self-consistent or self-consistent 
manner and (iv) the band structure method used. 

2.3 Review of Experimental Studies of the Fermi Surface 
of Vanadium 

The progress in the experimental study of the Fermi 
surface of vanadium is limited due to the difficulty in 
obtaining pure vanadium samples. The literature, therefore, 
contains somewhat fewer number of reports on experimental 
study as compared to the theoretical calculations. 

The earliest report on the experimental study was by 
Alekseevskii and Egorov [24] who carried out magnetoresis- 
tance measurements on V and observed an anisotropy in the 
FS of V, Later on Phillips [ll] used the technique of 
impulsive fields to study the FS of V in the tlOO } plane. 

This was followed by an extensive work by Parker and 
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Halloran [25] who used the magneto thermal oscillation 
technique to study the FS of V in the ^ 100 ^ and fllO) planes. 
These authors also compared their measurements with other 
experimental and theoretical results. 

General features of the FS topology for Y, as 
supported by different studies, are as shown in Fig. 2.1. 

We have also shown in Fig. 2.2 the cross section of the 
FS of V in the central planes {100 > and {110 > as calculated 
by Mattheiss based on (a) energy bands of tungsten and (b) 
energy bands of iron. 

The r -centered octahedral hole surface arising out 
of the second band is shown in Fig. 2.1 (a). The hole surface 
arising out of the third band is rather complicated and 
consists of multiply connected hole tubes along the <100 > 
direction (commonly known as the jungle-gym- arms) as well 
as large distorted half-ellipsoids centered at N (Fig. 2.l(h). 
The estimate of the dimensions of these i-T-centered half- 
ellipsoids is obtained from the measurements of Parker and 
Halloran [25] who reported the values of the semiaxes along 
N-P, N-r and N-H as 4.652, 4.405 and 3.649 nm*~^ respecti- 
vely. Although the measured areas of these half-ellipsoids 
[24,25] show a good agreement among themselves they appear 
to be higher than the theoretical values calculated by the 
APW method [ 13 ]. The experimental ratio of the semiaxes 
N-r /N-P, N-r /N-H for Y as reported by Parker and Halloran 
[25] are 0.946 and 1.21 respectively. In a similar way the 
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measured values of the effective mass (m*/m) determined from 
the de Haas-van Alphen experiments [ 25 ] are found to he 
higher than those calculated theoretically [13,18]. 

There is a major point of controversy about the FS 
topology of V. According to the model FS of V obtained by 
Mattheiss [7] by using the energy bands of Fe, the N- 
centered ellipsoids are connected to the jungle-gym arms by 
narrow necks extending along the r ~N direction. Although 
the experimental data by Parker and Halloran [24] admit the 
possibility of such necks, these authors were not able to 
confirm their existence. In terms of the band structure of 
V, the existence of such necks would require that the third 
band along the £ direction in V should not cross the Fermi 
level but should lie above it. Theoretical calculation of 
the band structure of V by different methods [13,15,18,20] 
report that the third band along the £ direction dips below 
the Fermi level by few- „ mRy. In view of this small 
'dipping' it is difficult to categorically confirm or deny 
the existence of these necks because the location of the 
energy bands near the FS is usually known within few 
milliRydbergs only. Alward et al. [19], for example, put 
an upper limit on the external diameter of these necks as 
2n/a ) if they exist. In Chapter 4 we shall discuss in 
detail the effect of these necks on the 2D-angular correla- 
tion of positron annihilation radiation. 
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2 . 4 Survey of the Compton Scattering and Positron 
Annihilation in V 

There have "been several experimental and theoretical 
studies of Compton scattering in vanadium. Weiss and 
DeMarco [26] used the e_/t 9 population ratio, obtained 
from the x-dif fraction paired reflections from Y, to 
calculate the directional Compton profiles of Y. These 
authors found a 8 % directional anisotropy in the Compton 
profiles j ^_qq (q = 0) and J^il (l = 0)* However the experi- 
mental anisotropy measured by other workers [27-29] was 
found to be small. 

A detailed band structure calculation of the direc- 
tional Compton profiles of V was carried out by Wakoh and 
Yamashita [16] using the KKR method and a state-dependent 
potential [153. Their theoretical results showed a fair 
agreement with the experimental profiles measured by 
Terasaki et al. [28]. In a subsequent work Wakoh et al. 

[30] repeated their calculations using an APW method and 
obtained results which were qualitatively same as before. 
Recently Itoh et al. [ 31 ] have measured the directional 
Compton profiles for V and VD Q gg single crystals. Their 
measured data for Y agree well with the theoretical calcu- 
lations of Wakoh et al. [ 30 ] - 

The angular correlation of positron annihilation 
radiation (ACPAR) in V has been measured by Shiotani et al. 
[32] and by Hasegawa et al. [33] using the long- si it 
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geometry. Experimental data .from both, these measurements 
showed a satisfactory agreement with the theoretical calcula- 
tions of Wakoh et al. [17]. 

A prominent feature about the experimentally observed 
directional Compton profiles [27-29] as well as the direc- 
tional ACPAR curves [32,33] from V is the existence of fine 
structure which can be understood in terms of the ES topology 
[17,30]. For example the Compton profiles as well as the 
ACPAR curves along the [lOO] direction show a dip or conca- 
vity at low momenta near p = 0. This structure has been 
ascribed to the closed hole pockets at point N and the 
multiply connected Jungle-gym arms of the third band FS. 
Measurements made on VD g gg indicate that introduction of 
deuterium fill up the hole bands of V and the Compton 
profiles [ 31 ] and the ACPAR curves [33] for VD^ g g show a 
disappearance of the fine structure (dip or concavity) in 
the [lOO] direction. 

These results clearly indicate that the electron 
momentum distribution in V is sensitive to the FS topology 
and that these effects can be probed by measuring the 
Compton profiles and ACPAR curves for single crystal speci- 
mens. It was pointed out that the 2D-ACPAR measurements 
can examine such effects in greater detail and can provide 
a wealth of information about the FS topology. In the present 
work we have, therefore, attempted to understand the effects of 
band structure and FS topology on the 2D- ACPAR curves in V 
and VD . 

<A- 
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The occupation number f(k,j) is identical to the Fermi-Dirac 
distribution and is zero for the empty state and is unity for 
the occupied state (at T = 0 K) of the electron. It may be 
pointed out that the positron wavefunction is calculated for 
the state ic, = 0 in its lowest energy band since it is therma- 

T 

lized before annihilation. 

The electron wavefunction .can be obtained by solving 
the Schrodinger equation* 

[_ V 2 + V($) - E(£)3 *£(r) = 0 (3.3) 

in the independent particle approximation using the standard 
methods. Although many such methods for solving the band 
structure problem are available [2-6], previous calculations 
of the 2D-ACPAR in metals have used the OPW method for A1 [7], 
the APW method for Ou [8] and the Hubbard-Mi jnarends approxi- 
mate method for Cu [ 9 ] and Fe [lO]. In practice the calcula- 
tion of the 21-ACPAR curve N(p ,p ) along a particular set of 
^Py»P 2 ^ direction requires computation of the TPMD, p(p), at 
about 500 p-points in a plane. 

This makes the computing time and computer memory a 
serious consideration. The method due to Hubbard and 
Mijnarends [ll,12] provides a fast but approximate scheme 
for such calculations and it has provided satisfactory 
description of the ACPAR data in transition metals [9,13*14]. 

We shall follow the mixed system of atomic units in which 
fi = 1, m = 1, e = 1, while the energies are expressed in 

Rydbergs. 
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These considerations along with the fact that these computer 
programs were available in our group guided us to use the 
Hubbard -Mi jnarends fast approximate scheme to calculate the 
TPMD in V. In the next section we shall outline the essential 
steps involved in this method. The detailed discussion of 
this method is given by Hubbard [ll], Hubbard and Mijnarends 
[l,12] and Harthoorn [l 53. 

3 . 2 The Fast Approximate Method due to Hubbard for the 
Calculation of Electron Band Structure ~ 

It was pointed out earlier that the electron band 
structure, wavefunction and the momentum distribution in a 
crystalline solid can be obtained by solving Schrodinger 
equation (3.3). A number of methods have been developed to 
achieve this [2-6]. The problem assumes a special character 
when we are interested in the electronic structure of transi- 
tion metals. The electron band structure of transition 
metals can be described as hybrids of nearly-free-electron 
(NFE) like bands with narrow d-bands of tight binding nature. 

A convenient approximate scheme using this approach and 
based on the Korringa [l 6 ] and Kohn and Ro stoker [ 17 ] and 
Ziman [6,18] method for the calculation of the band structure 
is given by Hubbard [ll]. In the Korringa-Kohn-Ro stoker (or 
KKR) method which is based on the Green' s function method, 
the crystal potential is usually assumed to be of the 
muffin-tin type i.e. spherically symmetric inside the 
inscribed sphere of the cell and constant within the 
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interstitial region. The band structure is determined using 
energy- depen dent phase shifts which arise from the scattering 
of the electrons from muffin- tin potential. 

In the KKR formulation by Ziman [l S ] the secular 
determinantal equation is given by the equation 

l [(k? - e) 6 , + r , ] B , = 0 (3.4) 

where = Ik^ } = |k + K n } , K n being the reciprocal lattice 

vectors and e is the Bloch energy eigenvalue for the electron 

in a state k. For the case where the potential term r is 

small, eq. (3.4) can be solved using only a relatively small 

number of plane waves. Lloyd [l 9 ] introduced the concept of 

r nn i , the matrix elements of the small pseudopotential which 

has the same scattering properties as the original lattice 

potential. He also introduced a pseudo-wavefunction _(r) 

ps 

which can be expanded in terms of the coefficients B n in the 
following manner 

* ps (?) = UA ) 1/2 l B n exp (i K n .?) (3.5a) 

throughout the atomic polyhedron and 

* ps A) = »g(5) (3.5b) 

in the region between the boundaries of the atomic poly- 
hedron (volume t) and its Inscribed sphere. 

Following Slater [ 20 ] , the matrix elements r nn » 
of the pseudopotential can be expressed as 
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r 

nn 


t 


where 


(*«A) I (2 2.+1) r \ T t (r t ) 5 t ( Vt )3 i (k n' r ! } V 00s tin’) 


(5-6) 


tan rj 4 

k r £ j £ (k r 4 ) {j £ (k r £ ) - (tan tj £ ) n £ (k r £ )> 


(3.7) 


In the above k = fe, r] £ is the phase shift of the & th partial 
wave scattering from the atomic potential, ® nn i ts the angle 
between k^ and k , while r £ can be chosen arbitrarily within 
the interval r = 0, r = r^ = the radius of the inscribed 
sphere of the atomic polyhedron. Other functions in (3.6) 
and (3.7) are as follows : P £ (cos 9) the Legendre poly- 
nomial, j £ and n £ the spherical Bessel and Neumann functions 
respectively. 

The quantity y £ and therefore r , will be small 
except the case when j £ (k r £ ) or [j £ (k r £ ) - (tan p £ ) n £ (k r )] 
is close to zero for some value of £ . To avoid this situation 
Pendry and Capart [2l] suggested that r £ can be chosen at 
the zero of n £ (kr £ ) provided that there is a zero in the 
interval (0,r^), otherwise r £ = r^ . The case when r £ is 
different from r^ will arise only when %- 0. 

It has been shown by Hubbard [ll ,22] that in the case 
of transition metals a d-resonance occurs when passes 
through tc/2, rising rapidly from 0 to a within a narrow 
energy range at low energies and hence there is no way for 
V 2 to prevent the appearance of a corresponding singularity in 
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the denominator of Y2* Following the rule by Pendry and 
Capart, will be small except near the resonance in tan rjg' 
Assuming that only one resonance occurs for & = X(x being 

O 

equal to 2 for transition metals) at k = e one can write 

0 

[ 11 ] 


V = Y ' 

‘£ 


6 £A W 

75 

e Q -k 


(3.8) 


where ¥ is the width of resonance and y^ is the residue which 
is always small. 

Ziman' s determinantal equation can now be written as 


IIC + y - £ | !h \ / B '' 


) 


\ h* W 


- o 


(3.9) 


where the matrix elements of the various submatrices are 
written as 


TC , 

nn’ 


= k* 6 , 

n nn' 


(3.10a) 


V nn' = IT l (2Jl+1) r ? <1 ^ r „) ^(k^.rj (cos 9^,) 


C ' I z n i' u c '"n 


£ £ 


nn' 


V = - 4 ” r x (|5 1/2 ■£ (Vl > T xm 


and 


a. 


m 


Z 3 x (k n r x ) T* m (ip B n 


(3.10b) 

(3.10c) 

(3-10d) 


with m = - x , — X + 1 ,X 

The column matrix and ^3 in the eq. (3-9) describe 
the resonant and non-resonant part of wavefunctions respectively 
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Following the procedure of Heine [ 23 ] Hubbard was able 
to reduce eq. (3*9) to an equation of lower order. In this 
procedure various vectors k n are divided into tv/o sets; (i) a 
'preferred' set P and (ii) a 'remainder' set R. The vectors 
in the set 'P ? are chosen so that for the vectors in set ' R ' 
not included in set 'P' the quantity (e - k ) never becomes 
small in the energy range of interest. As a result only a 
few of the smallest reciprocal lattice vectors will be 
included in the set ' P ' . 

This partitioning of the k^ results in further simpli- 
fication of the problem. The matrix elements of \\/pR and 
V^RP are a -*- so sm all because y^ is small. Using the approxi- 
mation = 0 and neglecting V^rr an comparison 

2 f R ) 

with (k^ e) we can rewrite eq. (3.9) in two parts as 



and 



where the matrix 

fcy 

A mm' " e o ^mm‘ 



elements of the square matrix 


(3.11b) 
are given 


(4m) 2 

T 


r 


2 

X 


¥ 


I (R) 

n 


.2 


Ox ( k n r x) 


n 
k 2 - 


Y Xm <*n> Y 


* fk ) 
xm' ^n^ 

( 3 . 12 ) 
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with the sum 


(R) 


n 


"being carried out over all vectors of the 


remainder set Huhhard [ll] has given the following 


prescription to calculate the above sum 


(R) 


n 


One writes 


„(R) 


(S) 


+ 


1 - cos ( 21c r, ) 
mm' r v o X 1 

r. t, L 2k r 2 1 


"X u k 


o x 


Si(2k Q r x )] 
... (3.13) 


Where 


(S) 


is a summation similar to the second term in eq. 


(3.12) over a set 'S’ of the vectors k n in the first few 

shells of the reciprocal lattice excluding those belonging 

x 

to the set *P’ , Si(x) denotes the sine integral / (sint)dt/t, 

o 

is the volume of the unit cell in the reciprocal space and 
N = (4Ttk^/3T v ) = total number of reciprocal lattice vectors 

OK 

(s) 

included in the set 'P' and the sum E . The determinantal 
equation (3.11a) is similar to the form of equation used in 
the model pseudo-potential scheme [23,24]. The physical 
significance of the various submatrices in eq. (3*lla) is as 
follows. The submatrix \jj describes the nearly-f ree-electron 
(NFE) part of the s-p conduction band while the submatrix ,;Al 
(5x5 in size) describes the d-bands. The off-diagonal matrices 
jb| and |«j* contain the hybridization between these bands. 

The energy dependence of the submatrices and ^//required 
further suitable approximation to remove the energy dependence 


of the matrix 


m 


"^(O - ejj 


Bp! 


1 t * — f ; 

\ V 


= o 


(3.14 


which is an equivalent form of (3.11a). The whole aim is 
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to reduce the problem to that of diagonalization. In Hubbard's 
scheme this is achieved by setting 






(3.15) 


such that jjVjj ^ and are energy (e) independent and by 

substituting 



to obtain 



as an equivalent f--rm of (3.11a). 

It can now be shown, by comparing eq. (3.11a) and eq. 
(3*16) that the actual problem to determine the band 
structure has now been reduced to secular form which can be 
solved by standard methods of matrix algebra. 

In the above only a bare outline of Hubbard's method 
is given. A detailed description of the method is given 
elsewhere [11,12,15]. 


3.3 Calculation of Electron Wavefunctions ; Hubbard- 
Mi jnarend s Method 


In the preceding section we have given a brief 
description of the fast approximation scheme due to Hubbard 
for the calculation of band structure of transition metals. 
This scheme, does not contain any adjustable parameter and 
is known to provide an overall accuracy of 0.005 - 0.01 Ry 
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[ 25 ] when compared with other methods like the APW method etc. 
In 1972, Hubbard and Mijnarends [12] extended this scheme to 
calculate the electronic wavefunctions in a form suitable 
for the calculation of matrix elements such as those in eq. 
(3*2). Using the concepts developed earlier these authors 
showed that the eigenvectors of eq. (3.4) contain the Fourier 
coefficients of a plane wave expansion of a pseudo-wavefunction. 

When a muffin- tin type of crystal potential is used, the 
electron wavefunction 4jt(r) can be expanded in two different 
forms in the two regions. The pseudo-wavefunction is diffe- 
rent from the true wavefunction in the region inside the 
inscribed sphere. In this region the true wavefunction can 
be expanded in terms of spherical harmonics Y^Cr) with L as 
the quantum number L=(£, m) together with those solutions 
R^Cr) of the radial Sehrbdinger equation which are regular 
at the origin. The electron wavefunction 4^(r) within the 
inscribed sphere is given by 


*£ (?) 


4x y . % 

( t ) 1 /' 2 1 1 MV 


X L Y L 



(3.17) 


In the interstitial region, the true electron wavefunction 
coincides with the pseudo-wavefunction and it can be 
expanded in terms of plane waves with B n as coefficients 

= (1A) 1/2 l B n exp[i(k+E n ).r ] 


*£ (r) 


(3.18) 
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The coefficients in eq. (3.17) can be determined by 
requiring continuity of ^ £ (r) at r = r . to obtain 



The next job is to normalize the wavefunction to 
unity in the atomic cell. Rigorous treatment shows that the 
wavefunction ^ (r) will be correctly normalised if the 
eigenvectors'^ in the secular equation (3.16) are correctly 
normalized. 

Although it is necessary to use all the jt-values 
up to infinity to achieve a perfect matching at r = in 
practice one limits i to fc usually say & „ = 3. This 
leads to a slight discontinuity at r = r.. 

3 • 4- Calculat ion of Positron Wavefunction 

In the independent particle model the crystal potential 
seen by the positron in a metal is same as that for the 
electrons except for the different sign. The facts that 
positron is thermal ized before its annihilation and that the 
average number of positrons in the sample at any given time 
(with the positron beam intensities commonly used in an 
experiment) is of the order of one further simplify the 
problem to that of determining the positron wavefunction in 
the lowest energy state k + = 0. 

Early calculations [26] of the positron wavefunction 
used the Wigner- Seitz approximation of a spherical cell. 
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In this method the positron wavefunction if» (r) = R + (r)/r is 
obtained by numerically integrating the radial equation 


d n. , (r) 


+ 


dr 2 


+ 2(E+V(r)) R (r) = 0 


( 3 . 20 ) 


and by choosing the eigenvalue to satisfy the boundary 
condition dif> + (r)/dr = 0 at r = r^g, the radius of the 
Wigner-Seitz cell. 

In view of the periodicity of lattice which is present 
in the positron wavefunction, it can be expanded in a 
Fourier series in terms of the reciprocal lattice vector [27] 


♦ .(?) = (1A) 1/2 l Si e 1K - ? 

K & 


( 3 . 21 ) 


The Fourier coefficients of the periodic 

potential V(r) may be found from the relation 

i(K-f ' ).r 




1 


Jv(*) 


at 


( 3 . 22 ) 


Now by substituting (3.21), (3.22) into the Schrodinger 
equation for positron, one obtains the secular equation as 

l {(K 2 -E) 6**, + Bi = 0 (3.23) 

i V j xi. Ja. y A. Ja. 

This equation is easily solvable after symmetrization. 

The Fourier coefficients are adjusted following Gould 
et al. [27] so that ! + (r = 0) = 0 subject to the normaliza- 
tion condition 

I , + , 2 

t 1 B £ I =1 (3.24) 

E 
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The computational details for the calculation of 
positron wavefunction in the case of bcc metallic vanadium 
is given in Chapter 4. The positron wavefunction along the 
three crystal directions, i.e. <100 > , <110 > , and <111 > can 
Toe calculated using the Fourier coefficients. The calcula- 
tion of positron wavefunction along the three symmetry 
directions offer the advantage of calculating the spherically 
averaged wavefunction R /r, by using the expression due to 
Houston [28]. 


R + (r) 


_ 1 _ 

35 


(10 


^100 ^ 


-f 16 




+ 9 


*111 


••••«• (3.25) 

This spherically averaged positron wavefunction was used in 
the calculation of the two-photon momentum distribution. 


3.5 Many -body Effects 

In the above discussion we have ignored the positron- 
electron (e’ r - e - ) and electron-electron (e~ - e”*) correlation 
effects and have considered the positron annihilation in the 
independent particle model. In reality there is a pile-up 
of electron cloud at tire positron due to the Coulomb attrac- 
tion thus raising the electron density at the positron site 
and enhancing the annihilation rate. Although these effects 
are significant for the positron lifetimes, they can be 
ignored for the angular correlations [l]. Nevertheless the 
enhancement of the TPMD arising out of the e + - 


e interaction 
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can be included in the following manner suggested by Kahana 
[29]. One calculates the TPMD P(p) on the basis of an 
independent particle model (IPM) and multiplies it by an 
enhancement, factor e(p) so that 

p(p) = e (p ) p(p) (3.26a) 

enh. 


where 


e(p) = a + b y + c y* 


and y = p/p ■ < 1 (3.26b) 

P 

and the constants a, b and c depend on the electron density 
[29]. Recently Mijnarends and Singru [9] have proposed an 
energy-dependent form of 'Kahana enhancement' which provides 
an excellent agreement between the theory and experiment for 
Cu. 

In the present work we have limited ourselves to the 
calculation of the T.PMD p(p) in the independent particle 
model only. Any correction to P(p) for the enhancement can 
be made by a suitable formula later on when a comparison 
between theory and experiment is desired. 

3 * 6 Calculation of the Two-photon Momentum Distribution 

As discussed in Sec. 3.1 it is necessary to evaluate 
the overlap integral A^(k,p) given by eq. (3.2). lo do so 
the atomic cell is split up in two regions ; (i) the inscribed 
sphere and (ii) the interstitial region, i.e. the remaining 
part of the 'Wigner-Seitz cell. For the integration in the 
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inscribed sphere, the equation (3.17) is substituted for 
electron wavefunction and spherically averaged positron wave- 
function given by eq. (3.25) is used, while the product term 
exp (-ip. r) is expanded in partial waves. The integration 
outside the inscribed sphere is performed using the electron 
wavefunction and positron wavefunction expanded in plane waves 
according to (3.18) and (3.21) respe ctiveljr . The integration 
of this double series over the full Wigner-Sei tz cell then 
yields a double sum of 6-function. The integral of the plane 
wave expansion over the inscribed sphere must be subtracted 
from the last integral performed over both the regions. The 
final expression is 


-V -V 

A j(k,p) 


l 

n,n' 


B 


n 


B 


+ 

n’ 


[ 6(k +K 
L v n n 




+ 


4n 2 

172 


(v) 


l 

L 



r . 

i j jl (pr)R Ij; .(r)R + (r) r dr 
o 


(3.27) 


which formed the basis of our calculation. Some other 
computational details of our calculation are given in the 
next Chapter. 
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Chapter 4 


TWO -PHOTON MOMENTUM DISTRIBUTIONS AND 2D-ACPAR CURVES 
IN V: RESULTS* AND DISCUSSION 

4*1 Introduction 

In the present Chapter the results of our calculations 
for the electron hand structure, Fermi surface (FS) , positron 
wavefunction, two-photon momentum distribution (TPMD) and 
'two-dimensional angular correlation of positron annihilation 
radiation (2D-ACPAR) in vanadium will be presented and 
discussed. As explained in the previous chapter the starting 
point of our calculation is the choice of the muffin- tin 
crystal potential for vanadium. We have used the crystal 
potential prepared by Mijnarends [l] for the Special Project 
sponsored by the Special Commission on Electron Charge, Spin 
and Momentum Density of the International Union of Crystallo- 
graphy. This muffin-tin potential was constructed with the 

t 2 

help of neutral-free-atom charge densities for the 3d" 4s 
configuration obtained from a relativistic calculation. The 
Gaspar-Kohn-Sham exchange potential with a = 2/3 was used 
while the charge densities of atoms up to the 13th nearest 
neighbours were superimposed with the help of the Lowdin 

*The preliminary results of this work were reported at the 
Fifth International Conference on Positron Annihilation, 

Lake Yamanaka, Japan, April 8-11, 1979. The earlier work 
had used constant positron wavefunction. The present work, 
however, uses a realistic positron wavefunction in V. 
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Table 4.1. Numerical constants used for the crystal 
potential for metallic vanadium. 

P ar ame t er V al ue 

lattice constant* a = 5.74509 a.u. 

Radius*of the inscribed sphere r^ = 2.48770 a.u. 

Constant potential between 

the muffin- tin spheres V Q = ~ 1*7106 Ry 

Fermi energy with respect to the 

zero of the muffin- tin potential Ep = 0.746 Ry 

Mean d-band energy (the energy 
at which tan p , goes through 

resonance) = 0.839 Ry 

Width of the resonance W = 0,211 Ry 

*Taken from Mijnarends [l ] . 
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Table 4.2. 

The normalised muffin- tin potential [l] 

(r V(r)/2Z) at some selected values of r. 

SI. 

No. 

r 

in a.u. 

r V(r )/2Z 

1 

0.00000 

1.00000 

2 

0.00194 

0.99333 

3 

0.00389 

0.98659 

4 

0.00778 

0.97297 

5 

0.01167 

0.9594C 

6 

0.01945 

0.93283 

7 

0.02724 

0.90743 

8 

0.04280 

0.86063 

9 

0.05837 

0.81859 

10 

0.08950 

0.74545 

11 

0.12063 

0.68374 

12 

0.18289 

0.58433 

13 

0.24516 

0.50756 

14 

0.36968 

0.39711 

15 

0.49421 

0.31822 

16 

0.74326 

0.21103 

17 

0.99231 

0.14012 

18 

1.49042 

0.06129 

19 

1.98853 

0.02303 

20 

2.98474 

0.00132 
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alpha expansion. The potential was not male self-consistent. 
The parameters used in our calculation are given in Table 4.1 
while the numerical values of the normalized potential are 
given in Table 4.2. 

4 . 2 Band Structure and Band widths 

Using Hubbard' s approximation scheme described in the 
previous chapter energy band-values E(1c) for 91 S-points in 
the l/48th part of the Brill ouin zone (BZ) (Bigs. 4.1 and 
4.2) were calculated. Kaga [2] has pointed out the effect 
of the number (N ) of plane waves (the so-called 'preferred' 

ir 

set used to describe the nearly- free- electron part of the 
pseudopotential in Sec. 3.2) on the accuracy of Hubbard's 
fast approximation method [3]« Previous work [ 4 ] has shown 
that the effect of increasing from 7 to 19 was to shift 
the energies by about 0.005-0.03 Ry, the maximum shift 
occurring for the point We therefore used = 19 in 

the present calculations and these 19 'preferred' vectors are 
listed in Table 4.3. Terms upto l = 3 were included for 
describing the electron wavefunction inside the muffin- tin 
sphere. The band structure obtained by us along the important 
symmetry directions in V is shown in Pig. 4.3 while the mean 
d-band energy and the width of d-resonance is given in Table 
4.1. The density of states histogram, is shown in Pig. 4.4 
and the value of Permi energy (Ej>) obtained by us for V was 
Ej. = 0.746 Ry. The electron energy eigenvalues at some 





( 0 , 10 , 0 ) 



r (o ? o 7 o) 


F ( 3,7, 3 5 



A (3,3,3) 



?h| 91 k- points in the 1/ 4'- tn B Z of V for which 
K(k) values were calculated » 
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Tabl e 4 


—V 

3- The set of 19 preferred vectors (in units 
of 2n/a) used in our calculation. 


1 

( o, 

0, 

0) 

2 

( 0 j 

-1, 

1) 

3 

( 0, 

-1, 

-1) 

4 

( 0, 

-2, 

0) 

5 

( 1* 

-1, 

0) 

6 

(-1, 

-1, 

0) 

7 

(-1, 

0, 

-1) 

8 

(-1, 

1* 

0) 

9 

( 1, 

0, 

1) 

10 

( 1, 

0, 

-1) 

11 

(-1, 

0, 

1) 

12 

( 1, 

1, 

0) 

13 

( 0, 

1, 

1) 

14 

( 0, 

1, 

-1) 

15 

( 2 

0, 

0) 

16 

( ~2» 

0, 

0) 

17 

( o. 

2, 

0) 

18 

( 0, 

0, 

2) 

19 

( 0, 

0, 

-2) 
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. E { Rydberg } 
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Fig, 4»4r Density of states histogram for 7 , , The j 

..Fermi energy 2^ was 0.746 fiy. 
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Table 4.4. Electron eigenenergies of vanadium at some symmetry 
points determined by the present calculation. 


I .R. 

k (f) 

band 1 

band 2 

band 3 

band 4 

band 5 

band 6 

r 

0 0 0 

0.227 

0.795 

0.795 

0.795 

0.941 

0.941 

A 

o 

O 

0.447 

0.694 

0.822 

0.822 

0.915 

1.008 

H 

0 10 

0.505 

0.505 

1.054 

1.054 

1.054 

1.474 

Gr 

3 7 n 

TU IS 0 

0.540 

0.650 

0.766 

0.856 

1.075 

1.182 

N 

1 1 n 

7 2 0 

0.467 

0.671 

0.823 

0.944 

0.973 

1.086 

z 

3 3 n 

IS TO 0 

• 0.385 

0.710 

0.732 

0.928 

0.956 

0.965 

D 

113 

1 2 IS 

0.578 

0.649 

0.710 

0.966 

0.976 

1.285 

P 

111 

2 2 2 

0.652 

0,652 

0.652 

0.976 

0.976 

1.282 

A 

3 3 3 

10 10 10 

0.475 

0.672 

0.672 

0.955 

0.955 

1.092 

E 

3 3 3 

10 10 10 

0.625 

0.647 

0.647 

0.966 

0.966 

1.360 








Table 4.5. Comparison of energy band widths (in Rydberg) 
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important ic-points and. some important bandwidths are given 
in Tables 4. 4 and 4.5 respectively and the bandwidths are 
compared with the results of other calculations [5-ll] and 
experiments [12,13]. 

The bandwidths and band structure obtained by us show a 
satisfactory agreement with the results of Papaconstantopoulos 
et al. [5], the minor disagreement could be due to the diffe- 
rence in the potential, method of band structure used and the 
self-consistent character of the other calculation. In 
particular the second lowest band, which is interesting 

from the point of view of possible necks along the r N direc- 
tion (Sec. 2.3), dips below the Fermi level E = E p by about 
12 mRy as compared to 18 mRy as found elsewhere [14] . The 
agreement between our results for the bandwidths and those 
by Wakoh and Yamashita [6], however, is not so good although 
the band structure results look similar. 

4 . 3 Results for the Fermi Surface 

Using our band structure results and the value of 
E p = 0.746 Ry (Fig. 4.3) the Fermi surface (FS) of V was 
obtained as usual. The central cross sections of the 
FS in the (100 ) and (110 ) planes of V are shown in Fig. 4.5. 
Some other cross sections of the FS will oe presented later 
in this Chapter when we shall discuss the 2E-ACPAR curves 
in the p^. ^ 0 planes. 
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These results display the following familiar features 
of the FS of V : (i) Second hand hole surface centered at r f 

(ii) Third band distorted ellipsoidal surface centered at N 
(hereinafter called as NCE) , and (iii) third band multiply 
connected ' ' jungle-gym' 1 arms along <100> directions (here- 
inafter called as JG-A) . A rough estimate of the size and the 
shape of the NCE is given by the ratio of the semiaxes along 
NP, Nr and NH directions. Our results show that Nr/NP =0.88 
and Nr/NH = 1.16 as compared to the crude estimates 

Nr/NP = 0.95 and Nr/NH = 1.21 obtained experimentally by 
Parker and Halloran [l 5 ] using magneto thermal oscillation 
technique. The FS topology (Fig. 4.5) obtained by us agrees 
well with the results of Papaconstantopoulos et al. [5], who 
also find, in agreement with us, that there are no narrow 
necks connecting the NCE to the JG-A along the rN direction. 

4 . 4 Positron Wavefunction in V 

The positron wavefunction was calculated by expanding 
[16] it in terms of the Fourier series (Sec. 3.4) 



where s are the reciprocal ' lattice vectors. In the present 
calculation we used 8 shells of the reciprocal lattice vectors 
besides the vector (0,0,0). The positron wavefunction was 
thus expanded in terms of 141 plane waves. The same muffin- 
tin potential as in Table 4.2 but with a changed sign was 



Pos 
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H G N 



Pig. 4.7 j The rectangular grid of the k-pointa 
used i n^c al j^Jja ti ng the TPMJU, ? v\"p) f 
at the p(= k+iM ) points in the (lOO) 
and (1X0} planes. 
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used to calculate the coefficients Bj* following the procedure 
explained in Sec. 3*4. The lowest energy eigenvalue of the 
positron (at lowest point, k + = 0) was found to be 0.506 Ry 

with respect to the zero of the muffin-tin potential. The 
positron wavefunction, ^ + (r) along the three crystal direc- 
tions <100 > , <110 > and <111 > are shown in Big. 4.6. In the 
actual calculation of the TPMD , . p(p) , ( ec . 3.27) we used 
R + (r)/r, the spherical average of the positron wavefunction 
calculated inside the inscribed sphere with the help of the 
following expansion [17] 


R + (r)/r = [10 ^ 100 (r) + 16 ^-^(r) + 9 ^(r)] 

(4.2) 


4 . 5 Two-photon Momentum Distribution 

X y -y x 

The two-photon momentum distributions , o^(p) ,due to 
the band electrons in V was calculated at about 1200 p-points 
in the (010) and (110) planes using equations (3.1), (3.2) 
and (3.27) and a set of 141 K vectors. The mesh of the 

— y 

k-points used in these computations is shown in Fig. 4.7. 

The results for the TPMD, P^(p) , for the three symmetry 
directions <100> , <L10> and <111> are shown in Figs. 4.8 
to 4.10. These figures also show the partial contribution 
to the TPMD, P^(p), due to individual bands while the 


^Hereinafter we shall denote the TPMD due to the band 
electrons by the symbol p- b (^"). 
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relevant band structure is shown in the inset. The dotted 
curves show the value of the overlap integral, above the Fermi 
energy although thqy (being unoccupied) do not contribute to thi 
TPMD. 

As pointed out earlier in Sec. 3.1, the total TPMD, 

P ^(p) > due to the band electrons is built up from the contri- 
butions from the jth band through the summation 

occ 

V 5 > ■ ‘vtp) {4.3) 

1=1 


where ^ (p) is the partial TPMD due to the jth band 

electrons and the summation .is carried over the occupied 

bands. At a given point p = 3c + K. the value of Pt. H (p) 

o 

is determined by {A.(ic,p)| whose value, in turn, is deter- 

J 

mined by ^ . (r) and ^(r) (Sec. 1.4) as well as by the 
so-called symmetry selection rule. A group-theoretical 
treatment of the symmetry properties of the TPMD by Mijnarends 
[l8] and Harthoorn and Mijnarends [l9] has shown that the 
contribution to the TPMD from an individual band may possess 
nodes on the lines or planes of symmetry. This can be under- 
stood from the following expression 

-ip.? 


A 3 -(ic,p) = / e 


’l'J^(r) «l» + (r) dr 


(4.4) 


cell 

The behaviour of A ^(k,p) under the operations R of 
^ (k), the group of k, determines the symmetry properties of 
a ft- ^ -h^on « fun^tinn for the irreducible representation 
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^ o; ^ \fo ^ s "^e wavefunction . (?) and the term exp (-ip.?) 

should contain a part which will transform according to the 
same representation if the integral in eq. (4.4) has to be 
non-vanishing. This criterion is facilitated by the defini- 
tion of a projection operator 

9 5 = ( Vb) l X 3 (R)*P r (4.5) 

J R K 

where is the dimension of -r 5 , h denotes the order of 
( k ) f x^(R) describes the character of the operator R for 

this representation and P^ is the operator corresponding to 
R. When the projection operator, cp^ , is operated on a 
function the part belonging to ^ will be projected 

out of 0+. In case exp (-ip. r) = 0, the exponential term 

ir 

does not contain such a part and we shall have A^(k,p) = 0 
for this ff. For the value of A^(S,p) to be non-vanishing we 
should have <p ^ exp(-i^.r) o. This symmetry selection 
rule, therefore, plays an important role while determining 
the contribution of a particular band to the TPMD, particularly 
for a transition metal like V whose Fermi energy crosses the 
d bands. These effects are discussed by Harthoorn and 
Mi jnarends [l9] in some detail and we reproduce in Table 4.6 
some of their results because these are relevant for discu- 
ssing the behaviour of TPMD for V, It may be pointed out 
that the positron wavefunction does not affect the symmetry 
selection rule because the positron wavefunction is assumed 
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Table 4.6. Observable energy bands which give non-zero 
contribution to the TPMB, P-^p), at momenta 
for the bcc lattice [ 19 ]. 


p x Py p z Observable energy bands 
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To summarise, therefore, the partial contribution 
^ j(p) "by a jth hand to the TPMD is determined by the 
following factors: 


(i) The relation p = ic + , 

(ii) Occupied or unoccupied nature of the jth band at k, 

(iii) The electron wavefunction ^ . (r.) and the positron 

& 

wavefunction 4\(r), 

"T 

(iv) Whether the particular symmetry of the jth electron 

. ">■ *> . 

band allows it to contribute to A^(k,p; or not. 


These considerations have to be kept in mind while 
analysing the results of the TPMD in terms of the electron 
band structure and PS. In particular it has to be remembered 
that a PS discontinuity in the TPMD cannot always be observed 
in the first BZ (where it can be strong) because of the 
symmetry selection rule. In the higher zones the selection 
rule may allow observation of these continuities but with 
much weaker amplitude. 

The TPMD curve, P^Cp) , in V along the <100 > direction 
(Fig. 4.8) shows a smooth momentum dependence without any 
sharp discontinuity. This behaviour arises out of the fact 
that only the ^ bands are allowed to contribute along the 
<100 > direction to the TPMD by the symmetry selection rule 
(Table 4.6). Out of the two ^ bands (Pig. 4.8) only the 
lowest A band remains occupied throughout vi.e. from r to 
H along A ) and hence there are no PS discontinuities in the 
TPMD which is smooth. The boundaries of the PS sheet along 



the <100 > directions (i.e. the T-centered octahedral hole 
surface) arise out of the crossing of the band with the 
Fermi level and since A g land does not contribute to P^p), 
this 13 sheet is not able to produce a break (or discontinuity) 
in the P t,(p) curve of Fig. 4.8. 

The TPMD, P-^(p) , along the <110 > direction on the 
other hand shows sharp discontinuities (Fig. 4.9) arising 
out of the boundaries of the third band FS sheets (i.e. JGA 
and NOE described in Sec. 4. 3). Only the \ bands are allowed 
to contribute to the P^(p) along the <110 > direction and as 
a result the lowest two bands contribute to %(?)♦ 0ut 
of these two ^ bands the lower (or the ootsom— mos t) band 
is occupied throughout (i.e. from r to N along % ) while the 
higher (or third from the bottom) band remains occupied in a 
small central region midway between r and N where it dips below 
the Fermi level. These features are reflected in the P ^(p) 
curve (Fig. 4.9) which shows a spike-like structure around 
p^3.5 mrad. These FS discontinuities at p = 2.74 and 4.10 
mrad are repeated in the higher Brillouin zones giving rise 
to the so-called ’Umklapp’ images (or the so-called high 
momentum components) at p = 7*18 and 8.54 mrciU and p 14.03 
and 15.39 mrad respectively. Once again tne boundary of the 
r -centered octahedral hole surface is not able to cause a 
discontinuity in the P^(p) curve along the <110> direction 
because this boundary is caused by the V band whlch 1S 
not allowed to contribute to P b (£)* ae wlCltn of the Spike “ 
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like structure in the p b (p) curve (Pig. 4.9) is determined 
by the dimensions of the PS sheets JGA and NCE along the <110 > 
direction while its height is determined by the character of 
the "band below the Permi level. 

The TPMD, P-^(p), along the <111 > direction is shown 
in Pig. 4.10 and it shows a discontinuity around p* 11.10 
mrad. While going along the <111 > direction one follows the 
route r - A -p-F-H-F-P- A _ r etc. Between r and P only the A 1 
band is allowed to contribute (Fable 4.6) and the lowest A b 
band is occupied throughout while the higher band is 

unoccupied throughout. As a result the P- b (p) curve shows 
no discontinuities from F to P along A. On tne other hand, 
along the P direction (i.e. from P to H) only F^ hand is 
allowed to contribute and this F^ band (Pig. 4.10) crosses 
the Permi level midway between P and h. This crossing marks 
the boundary of the JGA between P and d ana it gives rise to 
the discontinuities in the P^Cp) curve (Fig. 4.10) at p *11.10 
mrad and again at p *16.55 mrad. It may be pointed out that 
the boundary of the JGA between r and P does not show up m 
the P-^(p) curve because it is caused by the band which 

is not allowed to contribute to p - b (p). 

The three curves for the TPMD (Pigs. 4.8-4.10) show 
features similar to the results for the MD (for the Compton 
scattering) calculated by Gandhi [ 20 ] . Although the PS 
discontinuities occur at the same points, the high momentum 

smaller than those in the EMD 


components in the TPMD are 






arbitrary ur*«ts 





Qfkiirari? units 
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for Compton scattering [ 20 ] . This behaviour is to he expected 

because of the effect of the positron wavef unction on the 
TPMD . 

In order to calculate the 2D-ACPAR curves we had to 
compute the line integral (Eq. 1.8) of the TPMD , p b (p) in 
the (010) and (110) planes along lines parallel to the p^-axis 
for a given orientation {p x ,Py,p 2 } . In order to provide a 
better intepretation of the calculated 2D-ACPAR curves we 
have plotted the contours (hereinafter called as isodenses) 
of the TPMD, P-^(p), in the two planes (010) and (110) 

(Pigs. 4.11 and 4.12). These curves bring out the highly 
anisotropic nature of the TPMD in V arising out of the 
complex PS topology and the d-character of the electron 
wavefunct ions at the PS. In particular the isodenses show a 
complicated nest-like structure around N and it arises out of 
the annular PS topology between the ITCE and JG-A (Pig. 4.5). 

The anisotropic nature of the TPMD is also brought 
out by the plots in Pigs. 4. 13 and 4.14 where + he TPMD, 

P b^x'Py = ^ > P z ) ™ Px P° r different values of 

P* ( viz. p_ = 0.0, 1.60, 3.19, 4.79, 6.39, 7.93, 9.58 and 
11.18 mrad) in the (010) and (IlO) planes respectively. These 
curves clearly show the effect of the PS discontinuties and 
A-(£,p 2 on the TPMD. 

4 . 6 Two-dimensional Angular Correlation 

The results of the TPMD, P^(p), due to the band 
electrons in V calculated for about 1200 p— points in the 
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Table 4.7. 

Crystal orientations {p 

z 

used for the calculation 

i.e. 2D-ACPAB. 

>p y > and 
of N(p y , 

momentum me sh 

,p_) curves 

Zx 

Orientation p . n 

2 * X 

Plane 


Mesh 5 


II 

o 



P z x p x 

A 

[001], [loo] 

(010) 


0.80x0.80 

B 

[ 001 ], [no] 

(no) 


0.80x1.13 

C 

[112], [ill] 

(no) 


1.96x1.38 

D 

[no], [no] 

(001) 


1.13x1.13 


5 In units of mrad (i.e. me x 10 ^ - 0.137 a.u.). 
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(010 ) and (110) planes on a rectangular momentum mesh (given 
in Table 4.7) were used to calculate the 2D-ACFAR curves 
from the following relation: 

N(p y ,p z ) = / p(p x ,p y ,p z ) dp x (4.6) 

The four crystal orientations used for such calculations are 
given in Table 4.7 and fig. 4.15 while the calculated 2D- 
ACPAR curves are presented in figs. 4.16 to 4.19- In these 
figures we have shown the derivatives of the K(p y ,p z ) curve 
in order to bring out the structure in these h(p ,p z ) curves. 
The bandwi.se (the lowest band is labelled 'a' while the 
label 'b' refers to lowest second or third band, depending 
on whichever has ( + ) symmetry) contribution to the N(p y »P z ) 
curves is also shown. While understanding the nature of the 
U (p jp ) curves it has to be remembered that the contribution 
by a particular band to N(p y ,p z ) is determined not only by 
the occupied nature of the band (i.e. fS topology) but also 
by the symmetry of the band (Sec. 4.5). I n case of the 
general ff-point in the (010) and (110) planes p^-aP) ~ ® 
unless the jth band has positive ( + ) symmetry. In view of 
this restriction we have plotted the fS topology of V m the 
extended zone scheme (figs. 4. 20 and 4.2l) in uhe planes 
(010) and (IlO). The thick curves in figs. 4.20 and 4.21 
are drawn through points at which band 2 and band 5 (counted 
in order of ascending energy such that the h conduction 
level is at the bottom of the first band) touch and exchange 






Fig* 4.21s San® 
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their positive ( + ) and negative (-) symmetry. The symbols 
2( + ), 2(-) etc. are used in figs. 4.20 ana 4.21 to indicate 
the symmetries of the second and third hand (since fourth 
and higher hands being unoccupied do not contribute) in the 
different regions of the (010) and (110) planes. It will be 
useful to refer to the isodenses (figs. 4.11 and 4.12) and 
the fS topology (figs. 4.20 and 4.21) while discussing the 
nature of 2D-ACPAR curves (figs. 4.16 to 4.19). 

The N(py,p z ) curve for the orientation A i p = 0, 
p z = [OOl] and p x = [lOO] shows prominent structure in the 
form of shoulders at p z ^2.4 mrad and p z ^ 5.C mrad (fig. 

4.16). It is seen that the contribution to K(p„,p_) by the 

y 2 

lowest (labelled 'a’) band falls smoothly with p,_ while that 

due to the higher (second or third band, depending on the ( + ) 

symmetry, labelled ' b ' ) band shows two peaks which give rise to 

the structure in the total N(p y ,p z ) curve, '"dien analysed in 

terms of the TPM I) curves (fig. 4.13'), the TP; ID isodenses 

(fig. 4.11) and the fS topology (fig. 4.20), the structure 

at d >2.4 mrad can be ascribed to the fall in P b (p) inside 

the NOB (i.e. the line of integration, parallel to the p x 

axis, passes through N-centered ellipsoid). It may be 

pointed out that for the low p-values there is no visible 

effect of the dimensions (or "boundaries) o± Ju-A on the 

i\T(p jp ) curve because in this region the hue co the 

y z a 

, third band is low (it varies as p for the d-states at 

small momenta). The shoulder in the K(p^,p s ; curve at 
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P 2 ^ mra d once again arises because the line of integra- 
tion (parallel to the p^-axis) has now moved out of the NCE. 
This behaviour also explains the two peaks in the partial 
contribution to the N(p ,p z ) curve by the band ! b’. The 
steady fall in the N(p ,p z ) curve after p z = S.O mrad is due 
to the general momentum-dependent decrease in the p^(p) as 
well as due to the discontinuities arising from the bounda- 
ries of the JGA around p_ = 5.6 mrad, and n = 6.4 mrad. To 

summarise, therefore, the N(p ,-p ) curve for the orientation 

y z 

A shows effects due to the MCE while the effects of the JGA 
are masked out because at higher (p > 2.5 mrad) momenta the 
boundaries of the JG-A lie in the 3(*~) symmetry region. 

The results of Fig. 4.17 show that the N(p y ,p z ) curve 
for the orientation B i p y = 0,, p z = [OOl] and p x = [110 J 
exhibits a concavity at p z v2.0 mrad. This concave structure 
is attributed to the NCE as well as the JGA and can be under- 
stood with the help of Figs. 4.12, 4.14 and 4.21. The 
position (i.e. the p z -value at which the N(p y ,p z ) is maximum) 
of the concavity is determined by the dimension of the NCE 
along NP . For example this position would move to a value of 
p z smaller than 2.0 mrad if the semiaxis of fne ellipsoid 
NCE (along NP) were smaller than the value obtained by us. 

The amount of concavity (i.e. N max - N(p z - = 0)), on the 
other hand, is determined by the dimensions along FN of 
the JGA as well as of the SCE. She hump-like structure 
(Pig. 4.17) at p *8-10 mrad arises frum the lowest band 
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’a’. It may be pointed out that a concavity is also observed 

in the Compton profile [21 ] and the long-slit ACPAR curve 

[22] from single crystal V for along the [lOO] direction 

z 

and these results are explained in terms of the ICE and the 
JG-A sheets of the PS. These experimental results involve, 
however, surface integration of P ^(p) and hence cannot be 
directly compared with our results which involve a line 
integral (4.6). 

In comparison to the orientations A and £, the I(py>p z ) 
curve for the orientation C : Py = 0, p z = [112] and p x = [ill], 
shown in Pig. 4. IS, is smooth and shows no structure. 

It is possible that some fine structure has been masked out 
by the relatively coarser (1.96 x 1.38 mrad"') momentum grid 
used for orientation C in our calculation. Another possible 
reason for this smooth behaviour is that the I® not 

symmetric about the praxis for this particular orientation 
C and any fine structure can get folded up while integrating 

P b (p) from p x = -“ to P x = +“ * 

The N(p ,p ) curve for the orientation D : p y = 0, 
p 2 = [llO] and p x = [HO] also show a concavity at p^l.O mrad 
(Pig. 4.19). This concavity arises out of the uCE as well as 
the JG-A as seen from Pig. 4.11 and 4.20. Eae position of 
the concavity is determined by the dimension of the PCS along 
the EH direction while the amount of the concavity is deter- 
mined by the dimensions along rl of both the JGA and the 
ECS. If the size of the ICS were to shrink along NH the 
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position of the concavity would move to p <1.0 mrad. The 

Cl 

depression in the partial curve ' b' (Pig. 4.19) in the region 
P z % 5. 5-7.0 mrad is caused by the iTCE. As a result, the 
curve (Pig. 4.19) for the total N(p y = 0, P z ) in the region 
?z ^ 7. 5-7 . 0 f all s off faster than in the case the NCE had 
not existed. 

As pointed out in Chapter 2 there is some controversy 
ah out the dimensions of the NCE and JGA along the TN (or z ) 
direction. The model PS of Y obtained by Mattheiss [23] as 
well as the measurements by Parker and Halloran [l 5 ] indicate 
a possibility that the NCE's are connected to the JOA’s by' 
narrow necks along rN direction. For such necks to exist 
the third band in V along the £ direction should not cross 
the Fermi level but lie above it. Our band structure result 
(Fig. 4.3) as well as those of Papaconstantopoulos et al.L5] 
and Wakoh and Yamashita [6] show the third band ( 2^) dipping 
below the Fermi level and thus preclude the existence of such 
necks. We have simulated these necks by lowering the Fermi 
energy below Ej, = 0.746 Ry by about 0.016 Ry (so that the 
third band is unoccupied) and have calculated the TPMD as 
well as the 2D-ACPAR curves N(p y = P z ) wi^li Ej, = 0.730 Ry 

for all the four crystal orientations A, B, C and I) (Table 
4.7). The results of these calculations are shown in Fig. 
4.2‘2 where we have plotted 

AN = N'(p y = 0, p z ) - N” (p y = 0, p z ) (4.7) 

(where N* (p = 0, p z ) and N» (p y = O r p z ) denote 
the values without neck and with neck respectively) 




Flit. 4.22*. 


Plot of tbs AN (eq. 4.7) curves to illustrate 
the effects, of the PS »eck oh, the 2D~ACPAh 

curves (see text). 
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for uie four orientations. The vertical "bars shown in Fig. 
4.22 represent 10 % of the N f (p y = 0, p z ) values at p z = 0 
for that particular orientation and thus give an idea about 
the relative magnitude of AiT. It is observed from Fig. 4.2 2 
that the presence or absence of the necks affect the 2D-ACPAR 
curves for each orientation significantly. Combining the 
results of Fig. 4.22 with those of Figs. 4.16-4.19 we see 
that the 2D-ACPAR curve for the orientation A will not show 
the structure at p z ^2.4 mrad and p z ^5.0 mrad if there are 
necks connecting the NCE's to the JG-A’s. Presence of these 
necks have no effect on the contribution to the 2D-A0PAR 
curve by the lowest ( + ) band (curve 'a' in Fig. 4.16; but 
it will affect the partial contribution (curve * b 1 in Fig. 
4.16; tty the higher ( + ) band. The effect of necks is to 
smooth out the curve ’ b T as a result of which the shoulders 
in the total ^(p, r >p,,) curve at p„ ^ 5.0 mrad will disappear. 
For the crystal orientations B and D, the important effect 
due to the presence of these necks would be that the concave 
structure around p = 0 (Figs. 4.17 and 4.19) would deepen, 

Zj 

the difference being maximum ( ^15 */* ) at p z = 0. The 
shoulder-like structure for the orientation D at P z * 7-9 
mrad would be smoothened out if the necks are present. The 
effects of the necks on the k(Py = 0, V%) curve for the 
orientation C, however, are not so sharp (Fig. 4. 22,.' . 

We thus find that the measurement of the 2D- AGP AR 
in V can throw light on the controversy of necks in on© 
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FS of V. It must be remembered, however, that in practice, 
the 2D-ACPAR curves are measured with a finite momentum 
resolution (of the measuring instrument). For example the 
effective momentum resolution of the multiple (64) ilal-detector 
system being used at the Brandeis University [ 24 ] corresponds 
to a two-dimensional resolution with nearly Gaussian cross 
sections of 0.6 mrad x 1.65 mrad FWHM (full-width-at-half- 
maximum) , Before comparing the calculated N(Py»P z ) curves 
with the experimental data it is, therefore, essential to 
fold the theoretical curves with the p - and p z ~resolution 
functions. The dimensions of the necks (connecting the 
NCB’s to the JGA' s) in V are expected to be very narrow, if 
they exist at all. For example, Alward et al. [llJ estimate 
that the dimensions of such necks would be^^(~) and would 
correspond to about 0*3 mrad for V* These numoers uierefore 
suggest that determination of the necks by 2D-A0PAR method 
is not so easy or straightforward. • Althougn we have 
discussed the case of V in the above discussion, tries e 
conclusions are equally applicable to Nb and Ta which have 

similar band structure and PS topology. 

^ q clear from the above discussion that the 2D-ACPAR 

% 

curves for V are very sensitive to the FS topology. Our 
calculations (described above) are limited to the four 
crystal orientations (Table 4.7) all in the p y = 0 plane. 

In view of the complicated FS of V one expects that 2D-A0PAR 

curves for V measured in the p y 4 ° P 1 ^ 13 ' as wel1 as for 
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Cross sections of the Fermi surface of V 1b di fl ©rent 
(p T p z ) * (110) planes, U) Plane jp, - 0* lb) 

p v «f2ra/10 t (o) Pv « Jl/2 1 H/10. la) * 3V2TH/10 
U) p„ , 4tf2fH/10 aM (f) py * 5I2TH/10, The- 
notation far tk* hatahed areas is sa&@ a® in Fig ,4,2 
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other crystal orientations) can reveal more detailed informa- 
tion about the PS topology of V. Although we could not calcu- 
late the TPMD and the 2D-ACPAR curves in the p =4 o planes 
because of the high cost of computer time involved, we have 
calculated the PS topology (from our band structure results) 
of V in p y =£ 0 planes parallel to the (010 ) and (110 ) planes. 
These results are shown in Figs. 4. 23-4. 25. While using 
these results (Figs. 4. 23-4. 25) for guessing the behaviour of 
the 2D-ACPAR curves lT(Py 7 ^ 0, p z ) two important points 
should be remembered (i) Although the hole FS sheets do not 
contribute to the TPMD, their boundaries (or dimensions) 
cause discontinuities in the TPMD which affect the 2D-ACPAR 
curves and (ii) In a general plane p 0 , all the occupied 
bands may contribute to the TPMD provided that the ( + )/(-) 
symmetry labels are not valid. 

The FS cross sections shown in Fig s. 4.23-4.25 emphasize 
the highly anisotropic nature of the FS of V. It is observed 
from Pig. 4.23 that the 2D-ACPAR curves H(p y 0, p z ) are 
expected to change sharply while going from p y = 0.1 rH 
(Pig. 4.23b) to p = 0.2 TH (Pig. 4.23c) and these changes 
would continue as p y is increased (Fig. 4.24). Such changes 
in the structure of the N(Py 7 ^ P 2 ) curves can provide 

valuable information about the various dimensions of the 
sheets of the PS. oimil&r argument holds for 
the results of Pig. 4.25* Ii ma y pointed out that some 
of the v = constant sections not shown in Pigs. 4.24 and 

rj 
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4.2 5 can be inferred from the symmetry considerations. We 
suggest that the measurements of the N(p 7 ^ 0 , p z ) should be 

planned while designing a 2D-ACPAR measurement on V. Results 
of such experimental measurements, when completed, can be 
interpreted by using the results presented in Rigs. 4.23 to 
4.25. 


4 . 7 Corrections for Core Annihilation and Enhancement 


The results for the 2D-ACPAR curves presented in the 

previous section describe the annihilations with the band 

electrons alone. At present there are no experimental data 

of the 2D-ACPAR curves in V reported in the literature. Ror 

a possible comparison between the theory (present work) and 

experiment at a future date our raw theoretical data (Sec. 

4.6) will have to be corrected. Apart from the correction 

for the finite angular resolution of the intrument (as 

discussed in Sec. 4 . 6 ) two other important corrections are 

necessary. Rirstly one should include the contribution to 

the 2D-ACPAR curves by the annihilations with the core 

electrons. In order to facilitate such a correction we have 

calculated the core TPMD, P c (p), due to the (Is ) 2 (2s ) 2 

( 2 n ) 6 ( 3 s ) 2 ( 3 p ) 6 electrons by the following formula [25] 

~ * rs 


p (p) = l (2*+l)l / 5i ( P r) R + (r! P nb r)r dr 

c n,& n 


(4.8) 


where jgpr) is the spherical Bessel fun< 


spherically averaged positron wavefnnotion and 


tion, R + (r)/r is the 
,(r)/r’ s are 



122 


’able 4.8. 


Two-photon momentum distribution (TP'MD), p (p), 

c 

of V contributed by the core [(is) 2 (2s) 2 (2p)' 
(3s) 2 (3p)^] electrons. 


SI. Ho. 


P 

( mrad ) 


P c (p) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


0.00 

0.0444 

0.73 

0.0447 

1.46 

0.0457 

2.19 

0.0472 

2.92 

0.0490 

3. 65 

0.0509 

4.38 

0.0527 

5.11 

0.0542 

5.84 

0.0550 

6-. 57 

0.0550 

7.30 

0.0542 

8.03 

0.0524 

3.76 

0.0497 

9.49 

0.0462 

10.22 

0.0421 

1C. 95 ■ ‘ 

0.0375 

11.68 

0.0326 

12.41 

0.0277 

13.14 

0.0229 

13.87 

0.0184 

14.60 

0.0144 


21 
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Table 4.8 (Continued) 

SI. No. 

. P 

(mrad) 

p c (p) 

22 

15.32 

0.0109 

23 

16.05 

0.0080 

24 

16.78 

0.0057 

25 

17.51 

0.0039 

26 

18.24 

0.0025 

27 

18.97 

0.0016 

28 

19.70 

0.0010 

29 

20.43 

0.0006 

30 

21.16 

0.0004 

31 

21 . 89 

0.0003 

32 

22.62 

0.0002 

33 

23.35 

0.0002 

34 

24.08 

0.0002 

35 

24.81 

0.0002 

36 

25.54 

0.0002 

37 

26.27 

0.0002 

38 

27.00 

0.0001 

39 

27.73 

0.0001 

40 

28.46 

0.0001 

41 

29.19 

0.0001 
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T abl e 4.9. 

Results for 

the 2D-ACPAR , N(p y = 

0, p z ) curves 


for V due to 

the isotropic core 

electron TPMD 


(Table 4.7) 

» c (p). 


SI .No. 


m N(p 

=t Py + P|) 



(mrad) 

1 . 


0.0 

0.0924 

n 

c. * 


0.5 

0.0908 



1.0 

0.0911 

4. 


1.5 

0.0911 

5. 


2.0 

0.0910 

6. 


2.5 

0.0908 

n 

i * 


3.0 

0.0906 

8 . 


3.5 

0.0901 

9. 


4.0 

0.0894 

10 . 


4.5 

0.0884 

11 . 


5.0 

0.0870 

12 . 


5.5 

0.0853 



6.0 

0.0832 

13. 


0.0807 

6.5 

14. 



15. 


7.0 

0.0778 



n c 

0.0745 

16 . 


i » J 


17. 


8.0 

0.0708 



8.5 

0.0657 

18. 



19. 


9.0 

0.0624 

20 . 


9.5 

0.0579 



10.0 

0.0532 

21 . 
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Table 4.9 (Continued) 

SI .No. 

, P . 

(mrad) 

N(p = Yp®+p|; 

22. 

10.5 

0.0485 

23- 

11.0 

0.0437 

24. 

11.5 

0.0390 

25. 

12.0 

0.0345 

26. 

12.5 

0.0301 

27. 

13.0 

0.0260 

28. 

13.5 

0.0222 

29. 

14.0 

0.0187 

30. 

14.5 

0.0155 

31. 

15.0 

0.0128 

32. 

15.5 

0.0103 

33. 

16.0 

0.0082 

34. 

16.5 

0.0065 

35. 

17.0 

0.0050 

36. 

17.5 

0.0038 

37. 

18.0 

0.0029 

38. 

18.5 

0.0022 

39. 

19.0 

0.0017 

40. 

19-5 

0.0013 

41. 

20.0 

0.0010 



L> i »'"«• a f f Sf r h : P bSTg«tr 

l«?'thgS : l*Bia«°o yA ’ i » , ^ So sh °» a( B0lid ourTe) 


by tho baa# el«ctr 
ahowa( solid carve) 
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the free atom orbital wavefunctions for V obtained following 
nerman and Skillman [26]. These values of the core TPMD 
(assumed isotropic) are given in Table 4.8 and it is about 
5 7 * of the TPMD due to the band electrons at p = 0. The 
contribution by the core electrons to the 2D-ACPAR curves is 
given in Table 4.9 and is also shown in Pig. 4.26* It is 
observed that the core electron contributions are about 
12 /. of the ni(Py = 0, p z ) curve due to band electrons at 
p 2 = 0 for the orientation A. 

It has been shown by Mader et al. [27] that in the 
case of Al the calculated core distributions have to be 
reduced by a factor of 2 with respect to the conduction 
electrons for obtaining a satisfactory agreement between 
theory and experiment. Such a core reduction factor has also 
been used by Mijnarends and Singru [28] for Gu. The physical 
argument in favour of the core reduction factor is that the 
core electrons, being tightly bound, show a smaller polari- 
zability in the presence of positron. The work of Mijnarends 
and Singru [28] has shown that the choice of the core reduc- 
tion factor depends on the type of enhancement correction 
(discussed below) and goodness— of— fit between experiment 
and theory. 

The second important correction arises because of 
the many-body effects. As pointed out In Sec. 3.5 the 
TPMD P^(p") , calculated for the conduction electrons within 
the independent particle approximation can be multiplied by 
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A 

*T 


8 


Calculation 
Wigner- Seitz 


oJ l.'fche 2D- ACPAB Gurv 
Model ' 


es_ 


Jo y Using 


Modifi ed 


In the case of Al, Mader et al. [ 27 J have calculated 
the 2D-ACPAR curves using a simpler modified Wigner-Seitz 
(MWS) model first suggested by Berko and Plaskett [25]. In 
this MWS model a free-electron sphere is centered on each 
reciprocal lattice site in the momentum space. The TPMD 
is then assumed to be constant within such a sphere and this 
value is taken to he W(K^) such that 

occ 

wkp I lAg f (4.11) 

j=l i. j 

where At _ depends on reciprocal lattice vector K. and- hand 
i,3 _ + + 1 

index j and it is equal to A^(k = 0, p = K^) because p =ic + 

The weights, for each sphere, were calculated 

from the results of full hand structure calculations described 
earlier in this Chapter. The geometry of the momentum space 
with the various Wigner-Seitz spheres (radii 5*92 mrad for 
V) centered at the different points for the BCG V are 
shown in Pig. 4.27. For those p-v alues which lie within the 
common volume of more than one sphere the weights ¥(K i ) have 
been added. The contribution from the (3d) core electrons 
was taken into account by calculations similar to those 
described in Sec. 4.7. The 2D-ACPAR carves, U(p y ,p 2 ) cal- 
culated by the above procedure (MVS model) are shown in 
Figs. 4.28 and 4.29. These results show that although the 
MWS model calculation is close to the full band structure 
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calculations for A1 [27], it differs significantly from the 
full "band structure calculation for 7. This behaviour is 
expected because the PS of V is non-spherical and highly 
anisotropic. The high momentum components in the MV/S curve 
are higher because of the absence of the hybridization effects 
which are taken into account by a full band structure calcula- 
tion in a natural manner. The results of Figs. 4.28 and 4.29 
illustrate the inadequacy of the MWS model to predict the 
fine- structure in the 2D-ACPAR curves and thus emphasize the 
need to perform full band structure calculation for a transi- 
tion metal like Y. 


4.9 Summary 

To summarize, the results presented in this Chapter 
for the TPMD and the 2D-ACPAR curves underline the importance 
of the effects of band structure and P3 topology in the case 
of Y. The discussion of these results point out the limita- 
ti ons of the orientations A, B, 0 and D for studying certain 
parts of the FS sheets. The FS cross sections in the P y ± 0 
planes suggest that measurements of H(p., + 0, P z ) curves 
are necessary to determine different dimansions of the FS. 
Our results indicate how the possible existence of the 
necks (connecting the KCS end the JSA) can be confirmed by 
the 2D-ACPAR measurements. The present theoretical calcula- 
tion provide a basis for comparing future experimental 
results with theory. Various correotionsrequired before 
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our theoretical results can be compared with experiment are 
also discussed. A comparison of our full band structure 
calculations with a simpler modified Wigner-Seitz model 
calculation is made and it shows that full band structure 
calculations are essential. It is suggested that a compre- 
hensive measurement of the 2D-ACPAR curves for the single 
crystals of V be undertaken by experimental groups. 
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Chapter 5 


RESULTS EOR THE 2D-ACPAR IN VANADIUM DEUTERIDE VD Q 72 
5.1 Introduction 

Recently the metal-hydrogen systems have become an 
interesting subject for theoretical and experimental research 
because of their technological applications as well as 
because of the fundamental nature of the problem [lj. 
Different theoretical and experimental techniques are being 
used to determine the electronic structure and other physical 
properties of the metal-hydrogen systems, ’when hydrogen is 
dissolved in metals it is found to occupy interstitial sites 
in the host lattice with an accompanying expansion in the 
volume. The metal-hydrogen system, therefore, shews changes 
in the electronic structure of the host metal which can 
arise from the influence of the hydrogen as well as from the 
structural changes in the metal lattice. Although tne 
experimental techniques involving the measurement of density 
of states or magnetic properties have been in vogue over 
the past many years [2], recent introduction of experimental 
techniques like Compton scattering [3,4] and positron annihi- 
lation [5] (which determine momentum distributions) nave 
provided new probes to examine different theories. 

Historically, there have been two main theoretical 
models which have been used to describe the electronic 
structure of metal hydrides [2]. She first model, called 
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M anionic model" arose out of the similarity in the proper- 
ties of alkali hydrides (e.g. LiH and NaH) and the alkali 
halide o (e.g. Lil , NaCl) both of which are characterised by 
the valence state of H . The main assumption in the ’anionic' 
model is that in a metal hydride hydrogen accepts an electron 
from the host metal and the metal behaves like M + and 
hyarogen like H . The behaviour of the transition metal 
hydride, on the other hand, suggests that in a metal hydride 
(e.g. v r H) hydrogen donates its electron to the host metal 
(without changing the band structure of the host metal) and 
assumes a 'protonic' character. This has given rise to a 
'protonic' model and according to it VH system should be 
similar to a substitutional alloy Y-Cr following a rigid hand 
appro ximation. Almost all of the earlier experimental data 
were being interpreted in terms of either the 'anionic' or 
'protonic' model. The ab initio band structure calculations 
of the electronic structure of the metal-hydrogen systems 
carried out in the last decade have, however, shown that the 
electronic states in a metal hydride are much different from 
those predicted hy the 'anionic' or 'protonic' model, A 
review of these hand structure methods is given by 
Switendick [2]. The results of the APW calculations in 
the transition metal hydrides show that in going xronx a 
metal to metal hydride there is little effeco on the 
electron states having predominantly d-character while 
states having s-like character shift to lower energies. 
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changing the hand structure considerably. 

It is known that vanadium can absorb large amount of 
hydrogen (or deuterium) which goes to the interstitial sites 
of the bcc host lattice. The electronic structure of the 
V-K and V-D system has been studied experimentally using 
techniques like electronic specific heat [6], magnetic 
susceptibility [7], nuclear magnetic resonance [8] etc. and 
the results were analysed in terms of the protonic rigid 
band model. Measurements by Fukai et al. [9] using soft 
X-ray emission spectroscopy indicated that addition of 
hydrogen (or deuterium) causes changes in the d bands of 
vanadium and formation of bonding states involving H(D)-ls 
and V-3d orbitals. These experimental results along with 
the theoretical calculations by Switendick [2] have shown 
that the electron states having si-like character at the H 
(or D) site are modified and the density of states for VH X 
(or YD, ) is different from that predicted by the 'protonic' 
model . 

Recently the electronic structure of vanadium 
hydrides and deuterides has been investigated by employing 
Compton scattering [5,4] ana positron annihilation techniques 
[5j. The idea in these studies is to study the electron 
wavefunetion through their momentum distributions, lasser 
and I eng el er [ 3 ] have measured the Compton profiles ox 
vanadium and palladium hydrides and have shown that the 
'anionic' or •protonic model is not successful in explaining 
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tneir results. The hand structure calculations which incor- 
porate metal-hydrogen bonding states (created below the Permi 
level by the introduction of hydrogen in the host lattice) 
nowever give a better agreement with the experimental 
results. Measurement of the directional Compton profiles 
from the single crystals of V and YD Q gg by Itoh et al. [4] 
have also shown that a fraction of the electrons from 
deuterium goes in filling up the hole states of 7 while 
the remaining electrons form bonding states between V and D. 

Recently Hasegawa et al« [ 5 ] have carried out long- 
slit ACPAR measurements on the single crystal of vanadium- 
deuterium alloys VDq VDq an ^ VDq 68* ^ iese authors 

observed prominent changes in the long-slit ACPAR curves of 
cure V and YD Q 6g for the [lOO], [lio] and [ill] directions. 
The characteristic finestructure in the long-slit or 1D-ACPAR 
curves for the pure single crystal of V (arising out of the 
complicated PS topology) was found to be absent in YI) g ^g g . 

The changes in the observed long— slit ACPAR curves can be 
attributed to a gradual filling of the ’ ' hole' bands of f 
with the dissolution of deuterium. 

The results of YD„ by Compton scattering [3,4] as well 
as by ACPAR studies [ 5 ] establish the fact that the momentum 
distribution in Y is altered by the addition of deuterium. 
One should therefore expect that a 2D-ACPAR study of VD x 
(or YH X ) system should provide us with a greater wealth 
of information about the band structure, PS topology and 
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momentum distributions in vanadium deuterid.es (hydrides). 
With this aim in mind we have calculated the 2I-ACPAR curves 
in VDq rj 2 and compared these with the results for pure V. 

The method followed in our calculation is outlined in the 
next section. 


5 . 2 Method of Calculation 


The phase diagram of the V-D system is less complica- 
ted than that for the V-H system (Figs. 5.1 and 5.2). The 
a-phase of the V-H system covers a limited composition range 
at room temperature. For VD X the a— phase reappears at room 
temperature for x > 0,65. The studies of.positron annihila- 
tion in the VH system by Mantl and Singru [lO] have shown that 

A 

there is trapping of positrons by lattice dexects which 
affects the positron parameters across the phase boundaries. 

In order to study the 21-ACPAR curves without the complica- 
tions of positron trapping one should study the VD X system 
rather than the VH X system so that even at room temperature 
the a-phase continues as long as x> 0.65. We have, therefore, 
chosen the system TO 0-72 for which the 2D-ACPAR curves are 

calculated. 

Perhaps the best way to calculate the IPMD m 


TO 0 would be to make an ab initio calculation employing 
an AIA method [ll] or the CPA-KR method [12] to obtain the 
electron wavefunotions. problem of finding the positron 


nnfflni i cated because the affinities 
wavefunction becomes more complicate 

„ v have to be known. The e -impurity 

of positron for V and D have 








Pig 
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scattering effects may also be important. Thus the exact 
determination of the TPMD in ^ by such ab initio method 
would pose serious difficulties. We have therefore avoided 
this approach and have used an approximate method based on 
protonic rigid band model outlined below. It is realised 
that our method is naive and is open to objections. Never- 
theless it provides a first stage of our understanding of the 
TPMD and 2D-ACPAR curves in the VD X system. We believe that 
our present calculations can be useful to design a 2D-ACPAR 
experiment in the VD X system. 

The method of our calculation is the one used for 
Cu-Zn alloy earlier [l?3 and is outlined below. - The basic 


assumption underlying our calculations has been taken from 
the work of Kanhere and Singru [ 14 ] who studied the theore- 
tical electron momentum distribution in a series of 3d 
transition metals (V, Or, Fe, Ni and Cu) . Their results 
indicated that the partial contribution to the 3MD by the jth 
band, that is the P^C?) curves along a particular <hW> 
direction do not change much in going from a metal to its 
next higher neighbour having the same crystal structure 


(i.e. in going from V to Or). The total M>, wMch ^ 
obtained by a sum over' the occupied bands 


P b 


0 c c ^ 

(p) = l c b,j (p) 


(5.1) 


however shows major 
energy (Ep) * Using 


changes due to the change in the Fermi 
this idea and assuming the 'protonic model 
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and. a rigid band approximation we proceeded by the following 
steps for the calculation of the 1PMD and the 2D-ACPAR curves 
in the VDq ^ system. 

i) Following the protonic model it was assumed that the 
deuterium in VDq ^ donates all its electrons to the 
conduction band of vanadium. 

ii) Assuming the rigid band model, the Fermi energy (Ed, ) 
in VDq was calculated from the density of states, 
n(E)dE, and Fermi energy (E-p) in pure V by using the 
relations 


5.0 = 

;p 

/ n(E) dE 

— CO 

(5.2) 

5.72 = 

TP ? 

p 

L n(E) dE 

(5.3) 


In other words the density of states, n(n)dE,was held 
rigid in going from V to The Value of E F 

determined by us was 0.800 Ry. 

iii) Following Kanhere and Singru [14 J the partial IPMD 

p b was considered rigid in going fr ° m Y t0 VI) 0.72' 

The total IPMD ^(p) for VD Q>72 was determined at the 

1200 p-points by summing eq. (5.3) up to the new E| =0.800 
1^0 2D-ACPAR curves were next calculated 
for the four orientations (Table 4.7) as before 
(Chapter 4)- 
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iv) In our method of calculation there is no difference 

between a hydride and deuteride and our results could 
he equally applicable to VHq^. We have quoted our 
results as for because experimentally VT) 0 ?2 

is more suited for measurements. 

5 • 3 Results and Discussions 

The results of our calculations for VD Q ?2 are shown 
in rig . 5.3 and 5.4. The theoretical 2D-ACPAR curves for 
VD 0.72 are compared with those for pure V through the 
difference curves 

4tI = K ™ 0 _ 72 - H V (Py-PP (5.4) 

in Fig. 5.5 where the vertical bars indicate the magnitude 
of 0.1 Ny (p = 0, p z = 0). It was pointed out earlier in 
Sec. 5.1 that the Compton scattering [3,4] as well as the 
long-slit ACPAR data [5] from single crystal specimens of 
VD have shown that the curves from VD are less anisotropic 
as compared to those from pure V. These authors found that 
the behaviour of VD X was closer to Or than to V thus sugges- 
ting that the electrons of deuterium fill the hole states 
of V. The FS topology of the VD 0<72 calculated in the 
framework of the rigid band approximation and the'pro tonic' 
model is shoivn in Fig. 5.6 and it shows how the hole surfaces 
of 7 are shru nk by the donation of electrons from deuterium. 
It is seen from Fig. 5.6 that the second band is now fully 











Fig. 


The difference curves AN v®q. 
2D-ACPA11 curves for VJa.72 and 
orientations (a) A, {*) u, ko } 
The scale of ^HJ- 0 saBie aS 
in Figs. 4.16-4.19. 


s AN (®q. 5.4 ) between the 
VBq, 72 and V for the - ■ ' 
(h) 3, (c) C and \,d) D. 
same as that for utpysrC* P z 
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occupied, while third and fourth "bands are partially occupied. 
This filling out of the second and third "band must cause 
changes in the TPMD contributed by the second and third 
bands with resulting changes in the 2D-ACPAR curves. 

Although the general nature of the N(Py = 0> P z ) 
curve for the orientation A for YDq ^ is similar to that 
for pure V, the difference curve, AN, shows significant 
structure in the momentum range p 2 ^2. 5-10.0 mrad. This 
structure can be understood with the help of fig. 5.6. For 
the orientation 3, the concavity in the low momentum region 
of the N(p = 0, p z ) found for pure V is absent in YD Q>72 . 

V 

Although the ( + )/(-) symmetry of the bands is not shown in 
Fig- 5.6, the considerations of the FS topology suggest that 
the concavity should disappear in VD Q _ 72 because the ll-centered 
ellipsoids (HOE) have shrunk in size. This point was discu- 
ssed in Sec. 4.6. Similar considerations when applied to 
the orientation 3 explain the disappearance of the concavity 
in VD 0 72 - m any case the difference curves 4l! (Fig. 5.5) 
show significant structure which should be measurable. 

In order to examine the possible nature of the 


2D-A0PAR curves in the p y + ° planes of VD 0-72 " , ' 2 haYe 

calculated the FS topology in (p x P 2 ) P arallel t0 thS 

(010) and (110) planes. Ihese results for V3 0 . 72 , snown rn 

Figs. 5.7 - 5.9, =an be compared with similar results for 

A r> % _ 4 . 25 . Comparison of these 
pure V shown in Digs. • 3 

n j _-u 0 4- -hhp b eh a,vi oui? of the A ^ 
two sets of figures reveal that the bcha 
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curves (eq. 5.4) can show more drastic changes in going from 
V t0 VD 0.72 when the 2D-ACPAR measurements are made in the 
P y ^ ° plane. It is now known from the ah initio calcula- 
tion of the band structure of transition metal hydrides [ 2 ] 
that the rigid band approximation and protonic model are too 
simplistic in their approach. However until such ab initio 
calculations of the TPMD are made, the present results of 
our simple calculation can provide a first step of under- 
standing the dependence of the 2D-ACPAR curves in VD Q .72 on 
the band structure and PS effects. Recently calculations of 
the momentum distribution of electrons in concentrated 
random alloys Cu 1-X Ni x [ 15 ] and Ag Q>2 J»d 0 Q [l 6 ] have been 
reported. However these calculations refer to the case of 
Compton scattering because they do not take the positron 
effects into account. These effects would consist of the 
positron wavefunction (which is not constant), positron 
affinity for the atoms A and B (which may be different), 
the enhancement and core reduction corrections. Talcing all 
these effects into account for calculating a realistic 2 D- 
ACPAR curve in YD V will take more effort. Till such time 
that such calculations are made, the results of the present 
work may he used to analyse experimental results. 

5 . 4 Summary 

The results presented in this Chapter have shown 
that the 2D-ACPAR curves show measurable changes in going 
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from pure V to VD X . These changes are closely related to the 
changes occuring in the hand structure and PS of V after 
introduction of deuterium (hydrogen). It is suggested that 
measurements of 2D-ACPAR curves in V and VD„ (ensuring a 
proper choice of the phase as discussed in Sec. 5.2) be 
undertaken to examine the changes AH (eq. 5.4) in these 
curves. In spite of the simplistic assumptions made in our 
calculations for VDq our results can he useful to interpret 
these experimental data when they become available. 
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